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Abstract 


This thesis is devoted to the theoretical analysis of phenomena occurring 
in nanosystems under nonequilibrium conditions. The objects of the research 
are two-terminal nanojunctions consisting of quantum dots with metallic and 
superconducting electrodes. The main purpose of the research is to analyse 
time-dependent observables and correlation functions describing the system, 
e.g. differential conductivity, occupancy of quantum dots, and effective 
electron pairing on the dots. The nonequilibrium conditions which can be 
applied to the systems are a sudden change in the potential difference between 
the electrodes, a quench in the energy level of quantum dots, a sudden change 
in the hybridization of the quantum dot with the superconducting electrode, 
and the level of quantum dots oscillating in time. Under the influence of the 
mentioned perturbations, the evolution of the system takes place, which is 
noticed in observables and can be proved experimentally. The analysis of 
the above-mentioned quantities reveals the presence of beats in the charge 
current, a change in the direction of the flowing current, as well as a transition 
between the singlet and doublet states, a type of electron dot configuration. 
In periodically driven systems, we observe the appearance of harmonic states, 
an additional electron transport channels that contribute directly to the 
differential conductivity averaged over the period of the driving. These 
theoretical studies are based on Green’s formalism extended by Floquet’s 


theory and the Heisenberg equation. 


Streszczenie 


Niniejsza praca skupia się na teoretycznej analizie zjawisk zachodzą- 
cych w nanoukładach w warunkach nierównowagowych. Obiektami badań 
są dwuterminalne nanozłącza skonstruowane na bazie kropek kwantowych, 
których elektrody mają charakter metaliczny i nadprzewodzący. Głównym 
celem badań jest analiza czasowo zależnych obserwabli oraz funkcji korela- 
cyjnych opisujących układ, m.in. przewodnictwa różniczkowego, obsadzenia 
kropek kwantowych, efektywnego parowania elektronowego na kropkach. 
Warunki nierównowagowe, którym możemy poddać badany układ to m.in.: 
nagła zmiana różnicy potencjałów pomiędzy elektrodami, skok poziomu en- 
ergii kropek kwantowych, zmiana hybrydyzacji kropki kwantowej z elektrodą 
nadprzewodzącą, a także oscylujący sinusoidalnie w czasie poziom kropek 
kwantowych. Pod wpływem wymienionych zaburzeń zachodzi ewolucja układu 
zauważana w obserwablach, których przebieg może zostać sprawdzony ekspery- 
mentalnie. Analiza wymienionych wielkości ujawnia występowanie dudnień 
w charakterystykach prądowych, zmianę kierunku płynącego prądu, a także 
przejścia pomiędzy stanem singletowym a dubletowym, typem konfiguracji 
elektronowej kropki. W periodycznie zaburzanych układach obserwujemy 
pojawienie się stanów harmonicznych, dodatkowych kanałów transportu elek- 
tronowego, które mają bezpośredni wkład do charakterystyk przewodnictwa 
różniczkowego uśrednionego po okresie zaburzenia. Badania teoretyczne opier- 
ają się na wykorzystaniu formalizmu funkcji Green'a rozszerzonego o teorię 


Floquet’a oraz rozwiązywaniu równania Heisenberg'a. 
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Chapter 1 
Motivation 


In recent decades, nanodevice engineering has become a very well-known 
and developed part of the industry. The method of manufacturing nanostruc- 
tures offers the possibility of constructing variety nanodevices, information 
storages, transistors and other components needed in the high-tech branches. 
The breakthrough in the easy fabrication of semiconductors and semicon- 
ducting heterostructures enables the experimental verification of theoretical 
predictions for new physical phenomena in condensed matter physics. It is 
amazing how experiment goes hand in hand with theoretical predictions and 
the design of nanostructures. This remarkable development of condensed 
matter physics is evident in a growing number of publications in this field, as 
well as by its many new specialistic subdisciplines. 

Following such trends in physics, this thesis addresses the ideas of the 
condensed matter physics and is focused on the processes occurring in the 
superconducting nanojunctions. Since the late 1950s, the nanojunctions con- 
sisting of the superconductors and quantum impurities have been extensively 
studied [1,2]. The preliminary research has been focused on the paramagnetic 
impurities and their influence on the superconducting order. It has been 
shown that the impurities have negligible impact on the isotropic supercon- 


ducting state in the bulk materials [3]. However, it has been proved that the 
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magnetic impurities strongly influence the superconducting order [4-6]. In 
nanojunctions, it is known that in the boundary area between the normal 
metal and superconducting material there occurs a mixing of those phases [7]. 

Unfortunately, the experimental examination of bulk materials gives infor- 
mation about their macroscopic properties. To observe microscopic physical 
processes we need to focus on the study of the quantum impurities. In the 
1980’s, the quantum size effect has been observed by the absorption spectrum 
of excitons in quantum boxes |8] or in colloidal nanocrystals [9,10]. Those 
quantum boxes are epitaxially-grown quasi-zero-dimensional objects made of 
semiconducting materials. In 1986, the evidence of energy quantization in 
those objects has been shown and they have been named "quantum dots" [11]. 
Not much later, it has been presented that quantum dots are characterised 
by discrete peaks in resonant tunnelling spectra, indicating an atomic density 
of states in zero dimension [12]. Since then, new fabrication methods of quan- 
tum dots had been developed, e.g.: the electron beam lithography of narrow 
quantum wells made of GaAs/AlGaAs [13] or InGaAs/InP [14], the Stranski- 
Krastanow method [15, 16], application of carbon nanotubes [17,18]. The 
easy and well understood fabrication process has contributed to the discovery 
of a wide range of quantum dot properties. This knowledge contributed to 
the development of laser based on the quantum dots [19], photo-detectors [20] 
and also the incoherent light sources [21]. One of the milestone experiments, 
performed by Bonadeo et al., was the coherent control, by phase-locked pi- 
cosecond pulses, of the excitation in a GaAs/AlGaAs single quantum dot on 
a time scale that is short compared to the time scale of quantum coherence 
loss [22]. This work has initiated the field of quantum computation based on 
the optical properties of quantum dots. It has been shown that quantum dot 
experiments can satisfy some of the criteria of DiVincenzo [23] conditions, 
necessary for constructing a quantum computer, e.g.: state preparation [24], 
state read out [25,26] and coherent manipulation [27, 28]. 


Furthermore, the use of quantum dots in quantum information can be 
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based on the spin manipulation. It has been proposed that in the super- 
conducting junction with a quantum point contact or the quantum dot, the 
Bogoliubov quasiparticles can be trapped in discrete Andreev levels [29, 30]. 
In contrast to the semiconductor quantum dots this type of device has been 
called Andreev quantum dot. This proposal of Andreev qubit shows that the 
experimental manipulation and measurement of a state defined by spin 1/2 is 
possible, which provides the opportunity to quantum computing applications. 
Moreover, another theoretical concept of qubit proposed by Nazarov [31] 
is based on a superconducting circuit and quantum dot in the form of the 
Josephson junction. In systems with the spin-orbit interaction, the two-level 
system is realized by the spin state of the electron on the quantum dot and 
its polarization can be tuned by the superconducting phase difference across 
the junction. Based on those ideas, various families of qubit devices have 
been proposed and experimentally examinated [32]: the charge qubits in a 
transmission line resonator - transmons [33,34], the gate tunable supercon- 
ducting qubits - gatemons [35-37] or combination of transmon qubits with 
qubits based on Majorana zero modes [38-42]. This shows that quantum dot 
devices have a wide range of applications in the current and future quantum 
computing. 

Besides the single quantum dot devices, the double quantum systems have 
been experimentally studied by the tunnelling spectroscopy using InAs [43—46], 
InSb [47], Ge/Si [48] and carbon nanotubes [49,50] as well as by the scanning 
tunneling microscopy applied to molecules deposited on superconductors 
[51-55]. It has been shown that in such types of systems it is possible to 
manipulate and change the state of a single electron [56-58]. It had been 
proved that for the electrically isolated double quantum dot the eigenstates 
of the system are characterized by surprisingly long coherence time ~ 220ns. 
For quantum computation purpose the double quantum dot systems can be 
regarded as a kind of qubit [59]. 


Furthermore, superconducting circuits based on double quantum dot 
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systems are rich in the properties of in-gap bound states |47,60—62], originating 
from the ground state configuration with even or odd parity. Due to the 
system parameters, the parity change has been experimentally detected by 
discontinuities of the Josephson current in S-DQD-S junctions [45,46,63] or 
the subgap Andreev current in N-DQD-S junctions [47, 63,64]. In double 
quantum dot systems coupled to superconducting lead the spin of quantum 
dots could create a triplet blockade [63]. This spin locking counteracts 
the superconducting proximity effect, what has been reported in 5-DQD- 
S [44] and N-DQD-S [64] nanostructures. Quantum dot heterostructures 
can be studied in terms of time-dependent processes that occur under the 
influence of step-like pulses or photon exchange with the environment [65-67], 
leading to the post-quench evolution and photo-assisted tunnelling processes. 
The time-dependent experimental measurements for the Josephson junction 
placed in microwave cavity exhibits discretization of super-current and higher- 
order effects such as multiple Andreev reflections and Shapiro-like steps in 
conductance characteristics [68-70]. 

This thesis is devoted to theoretical studies of time-dependent processes, 
occurring in systems with single and double quantum dots placed between 
the normal metal and superconducting leads. The main objective is to study 
the evolution of systems induced by the temporary and long-term disturbance 
of its properties such as: rapid change of source-drain voltage or energy 
level position of quantum dot, sudden coupling to superconducting lead and 
periodic driving of the quantum dot energy levels. Depending on the type 
of perturbation, the systems exhibit different ways of reaching equilibrium. 
We investigate the relaxation times and associated phenomena manifested in 
post-quench evolution, evident in the time-dependent current characteristics. 
For this purpose we made use of: Green’s function method, Floquet formalism 


and Heisenberg’s equation of motion technique. 


Chapter 2 
Introduction to superconductivity 


The superconductors, as their name suggests, are characterised by perfect 
conductivity, which can be observed in specific materials under some physical 
conditions. Superconductivity, as physical phenomenon has been discovered 
more than one century ago [71]. The presence of perfect conductivity is one 
of the essential characteristics of superconductors, playing an important role 
from an economic point of view. Superconductors are also characterised by 
ideal diamagnetism [72]. In the presence of a magnetic field, its lines are 
expelled from the bulk of superconductor - the Meissner effect. For corre- 
spondingly higher values of the external magnetic field the superconducting 
effect disappear, and the field penetrates the volume of the material. The 
value of the magnetic field at which the superconducting effect is suppressed 
is called the critical point Be. It happens for the pure superconductors, like 
27 metals (e.g. Al, Pa, Th, Pb) which are called Type I superconductors, 
having the critical temperature To < 7.2K, below which the superconducting 
effect appears. Type II superconductors are made from alloys (e.g. NbTi, 
PbMoS, NbN), they are colloquially known as dirty superconductors. They 
maintain the superconducting state to higher temperatures (tens of Kelvin) 
and magnetic fields, than I type superconductors. The second type supercon- 


ductors can exist in a mixed state of normal and superconducting regions, 
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called a vortex state. This state predicted by Abrikosov is characterized by 
vortices of superconducting currents surrounding cores of materials in normal 
state. 

In this chapter, we will provide a short theoretical description of the 
superconductors and the processes occurring at the interface of different 
materials. The following sections (Sec. 2.1, 2.2) are based on the popular 


books [73-75]. 


2.1 Electron bound states of bulk materials 


The breakthrough step in the formulation of the microscopic theory of 
conductivity had been made by Cooper [76]. In 1956, Cooper showed that 
two electrons, interacting on the Fermi sphere, can create a bound state, 
called later a quasiparticle or Cooper pair. This state appears in the presence 
of the weak attractive interaction. One of Cooper’s assumptions was that 
pairing between electrons occurs in the singlet state X = wal N-I), 
which guarantees the antisymmetry of wave-function. In theory, the pair of 
electrons can be described by the two particle wave-function in the real space 
coordinates r; 


w(r1,12) = $ą(r1 — reje Ty, (2.1) 


where R = (rı + r2)/2 is the center of mass, $, is wave-function in the center 
of mass and q is momentum of the center of the mass of the pair with respect 
to the Fermi sea. One of the Cooper's assumptions is that the center of mass 
of the pair lies in the center of the Fermi sphere q = 0. Making a Fourier 


expansion of $,, the two particle wave-function takes form 
w(ri,r2) = ) g(k)ET te hy, (2.2) 
k 


where g(k) is a Fourier transform of the center-of-mass wave-function. The 
wave-function of the pair of electrons consists of the plane waves with equal 


and opposite momenta. 
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Solving the Schrödinger equation for the introduced wave-function (11, ra) 


of the Hamiltonian 
R h2 
HPS (WPN Vrs); (2.3) 


where V (71,72) is the interaction term, yields 


3 Fa + vo |o jett = (+25 +2 


k! 


Ny 2,y(k je (2.4) 


Here, the translational invariance has been assumed V(r1,r2) = V (p) and e 
is a total energy of two considered electrons. Multiplying both sides of (2.4) 


by ett? and integrating over d3p, we obtain 
2 
= f Dee rg ennratys | Dv era) dp= 
k! 
ace , 
(= i) p> D Weg(ki dp. (2.5) 


It is important to notice that the above equation has a non-zero solution only 


ifk=k’: 


7 d petre = Lox kr. (2.6) 
Then (2.5) takes form 
h2 
zk gk) + X` g(k View = ( + ŻEp)g(k), (2.7) 
k! 


where Vin = 73 f V (p) edp. 

The Cooper theory assumes that the attraction potential takes a constant 
value up to a cut-off imposed by the Debye frequency wp which is order of 
kp@p, where kg is Boltzmann constant and 0p Debye temperature. This 
means that the coefficient Vy, fulfills the following relations: 


-V/IB if |gx| Ew] < Aw, 
Vek = (2.8) 


0 otherwise, 
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where & = LC Substituting this simple form of potential Vi, into (2.8) 


2m 


we obtain 
V 1 
(-2& + e) g(k) => 3 2,9(k ). (2.9) 


Dividing both sides of (2.9) by the factor (-24;, +e) and performing summa- 
tion over k, yields 

oł) = 5 > g ak) (2.10) 
The right-hand side of (2.10) contains term Èx g(k’), which is generally a 
constant value, equal to the left side of the equation. This term can be 
omitted, which results in the self-consistent condition 


V 1 
Fis zl 


(2.11) 


Converting summation to integration from zero up to the Debye frequency, 
yields 
RER GE zy. (2.12) 
0 2E-E 
where v(€) is the density of states for electrons (which practically takes 
constant value over the range of the integration). Then, performing the 
integration and assuming that |e] « hwp, we finally obtain the energy electron 


pair 
E = —2hwpe HOW, (2.13) 


A negative sign indicates that a bound state is formed from attracting electrons. 
In other words, the energy of the bound electrons is lower than that of the 


free particles. 


2.2 Bardeen-Cooper-Schrieffer theory 


This scenario of the bound state formed of two electrons was fundamental 


to the macroscopic theory of superconductivity. Theoretical explanation of 
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the superconducting state was presented by a group of scientists consisting 
of J. Bardeen, L. Cooper and J. R. Schrieffer (BCS) and was awarded the 
Nobel Prize in 1972. The BCS theory is based on the Cooper model, the fact 
that the attractive interaction can bind two electrons with opposite spins 
and momentum vectors (k (, -k |) into one pair on Fermi sphere [77]. The 
Cooper result (see Sec. 2.1) implies that many-body system can reveal the 
energy gap in their excitation spectrum. For the many-body case, the Cooper 
result suggests that the ground-state wave-function Wy could be constructed 


of the paired electrons 


Wy = 9(1,2)Y(3,4)9(5,6)..9(N-1,N), (2.14) 


where N denotes the number of particle in the system. The pair wave-function 
W(1, 2) = Yy g(k)e*"ie"ikrzy is multiplied by individual coefficient g(k). For a 
realistic case this description is useless, because ~ 10? individual coefficients 
need to be calculated [74]. 

In the BCS theory and the second quantization representation, the ground 


state of the many-body system was proposed as collection of Cooper pairs 
|Y gcs) = a ky (Uk + vet ely )10), (2.15) 


where |0) denotes the vacuum and ug, vz, are coefficients whose square are 
the probability of empty (|u,|?) and occupied (|v,|?) state with wave-vector k 
respectively. In order to satisfy the normalization condition (Wpgos|WBcs) = 1, 


the probabilities must obey the following relation 
Juz |? + lvgl? =]. (2.16) 


The wave-function Wy defined in (2.14) describes the system consisting of 
N/2 electron pairs, while Wpos defined in (2.15) is a superposition of pairs 
containing the even number of electrons, without fixed number of particles. 


Between those two wave-functions a following relation can be written 


Wpos =) AwWn, (2.17) 
N 
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with the normalization condition È} yAn = 1. 


The many-body system with translational invariance can be described by 


Hamiltonian 
H=>&¢c 2 + YY, et el lać (2.18) 
m k-ko ko kk! ,oo'„qCk+q,o Ck'-q,o' Ck'o'Cko G 
ko kk',o0',q 
242 ę . . 
where ék = nk — €p denotes the electron energies with respect to the Fermi 


level. The second term of Hamiltonian in (2.18) describes electron-electron 
interactions. 

The binding energy of the Cooper pair depends on the sum of the total 
angular momentum of the electrons K = k,+k2, where kı and ko are their wave- 
vectors. It reaches a maximum value for K =0 and decreases dramatically for 
other values, implying that the electrons have opposite wave-vectors kı = —ko. 
The origin of the attractive force originate from a phonon exchange, giving 


rise to an effective interaction. Considering the negatively charged electrons 


is k+q 
Z 
q 


Figure 2.1: Diagrammatic representation of the exchange of a phonon q in the 


Eo 
I 


effective electron-electron interaction. 


moving through the lattice of the positively charged lattice ions. The phonon 
can be exchanged as schematically presented in Fig. 2.1. In this process, the 
maximum possible phonon energy is truncated by the Debye frequency wp 
(the crystal lattice is not able to transfer higher energies) and furthermore 
the angular momentum should be conserved k + k' = const. 

The Hamiltonian in (2.18) can be written in a simpler, reduced form. 


We can restrict the consideration to that part which contributes to the 
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superconducting effect, i.e., the interactive part taking into account only 
electrons with opposite spins and momenta. We can replace k' > -k and 


k+q—k', which yields 
H= Do + > | PRACZE (2.19) 
ko kk! 


To determine particular the coefficients w, and vy, we can use the variational 
principle of quantum mechanics to minimize the expectation value of the 
Hamiltonian (2.19) 
ó(WBesl Exel êl, + X Vine! ć-wyćwkbBos) = 0. (2.20) 
ko kk! 


The expectation value of the first and the second terms presented in (2.20) 


are given by 


(Vecs X &êt êl lbeos) = 2 Y Gzlvk|? (2.21) 
ko k 
and 
(wacs| X Verê êt p êr Cua laes) = )) Vek VZ UpUp Up (2.22) 
kh! kk! 


For convenience it can be assumed that the coefficient w, and vk are real 
numbers 


Up = cos and vy=sinQy, (2.23) 


which gives us the variational minimization condition with respect to the 


parameter 6, in form 


1 
za ( Y & sin” 6, + — D> sin 28, sin 20 =0. (2.24) 
005 k 4 ie 
As a result of minimization we obtain 
1 
4, tan 20; = =, Y Vir Sin 28w. (2.25) 
k! 


At this step, let us define an auxiliary parameter of the energy dimension 


1 
A = — 5 Vir URUK! = 5 > sin 20K", (2.26) 
k! k' 
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which can be called a gap function or order parameter. Substituting (2.26) 
into (2.25) we obtain 


A 
sin 20, = 2ugv, = —, (2.27) 
Ex 


where Ep = \/€? + Az. For completeness, we can also calculate 
cos 20; = uz — uz = Se (2.28) 


Having defined the relation between v and w, in (2.16) and their combinations 


in (2.27), (2.28), we can obtain the coherence factors in formal forms 


2 
v=5(1-=), (2.29) 


where the excitation energy Ep = yE? + Az. Substituting (2.27) into (2.26) 


yields a non-linear gap equation 


Ap 
An = >) Ver £ 


K IJE FA; 


The BCS theory assumes that the attraction potential takes a constant value 


(2.30) 


up to a cut-off imposed by the Debye frequency wp. This means that the 


coefficient Vx fulfils the following relation 


-V if + |Sk” < hw , 
Gene [ór] [Eve D (2.31) 


0 otherwise, 
which implies 


A for < ħwp, 
Az xl < hp (2.32) 
0 for RA > hwp. 


We can consider the following extreme situations: 


e |f| > A the energy of the particles is greater than the energy of the 


superconducting gap and they behave like free electrons, 
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e || «<A the electrons experience the pairing energy gap. 


Assuming momentum-independent potential (justified in isotropic systems) 
parameter A; can be cancelled on both side of the equation. Integrating 


(2.30) over energy € in the range +hwp yields 


1=v(0)V Vee ee = v(0)V sinh“! EBI (2.33) 
OT 
_ ney ae 
A = hwpsinh (or) ‘ (2.34) 


where v(0) is the normal density of states at the Fermi level. In the weak- 


coupling limit Vv(0) « 1 the gap function is given by the expression 
A 2 2ħwpe™ t! V”, (2.35) 


In the Hartree-Fock-Bogoliubov approximation, the interaction term in 


(2.19) can be linearized 
-V > êl êl pêkên = - a ENIA A Cnt) 
sd AC >, ĉr êr FV Pêl) 2 (ênir), (2.36) 
where V ¥°,(C_xyCxp) = A. It yields the reduced BCS Hamiltonian 


H= &é êl- (a Dalel + he] | (2.37) 
ko k 


Near the critical temperature the energy gap can be approximated by 


following relation: 
A(T) x A(04/1-— (2.38) 


where A(0) x 1.74kg7¢ is the energy gap for T=0. In the bulk supercon- 


ductor, the spatial range of the Cooper pairs is called the coherence length 
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Es and is strictly connected with the energy gap A. The coherence length 


ranges from tens to hundreds of nanometres and is defined as follows 


hur 


Es = TA’ 


(2.39) 


where vp denotes Fermi velocity. The temperature dependence of the co- 
herence length depends of the mean free path for the electrons le. In the 
homogenous materials és « le the coherence length has form 


T, 
T.-T 


£(T) =0.74€5 (2.40) 


For doped materials £s > le the coherence length obeys the following equation 


€(T) = 0.85,/Esle TT (2.41) 

In contrast to the phonon superconductors, we can mention other types 

of superconducting phases. For example, the superconducting state has been 
found in the heavy fermion compound UGeg [78], which is qualitatively similar 
to the classic d-electron ferromagnets. In the ferromagnetic superconduc- 
tor the quasiparticles responsible for the ferromagnetism form the Cooper 
pairs, for which the triplet pairing appears (k 1,-k f) [79]. Moreover, the 
superconducting effect can exist in form of FFLO phase (Fulde-Ferrell-Larkin- 
Ovchinnikov) [80-82]. In the presence of a magnetic field, we can find a pair 
of electrons (k t,q- |) in the bulk of the metal, which have nonzero total 
momentum K + 0. The external magnetic field induces a Zeeman energy 
splitting, the Fermi spheres differ for individual spins, resulting in the pair- 
ing of electrons into singlet pairs with total momentum q. Furthermore, it 
is important to mention about the high-temperature superconductors, the 
class of the materials which can not be simply described by the BCS the- 
ory. High-temperature superconductors exhibit precursor pairing effects at 
higher temperatures than T. Lately, it has been possible to obtain the first 
superconducting state at room temperature, discovered in the mixture of 


H>S+CH/ at high pressure 267 GPa [83]. Experimental results indicate that 
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the high-temperature superconductors involve singlet Cooper pairing with a 
d-wave symmetry [84], unlike in the BCS theory where the superconducting 


energy gap has a s-wave symmetry. 


2.3 The proximity effect and Andreev reflection 


At the metal-superconductor boundary the properties of both materials are 
mixed nearby the junction region [85-88]. The superconducting order enters 
the metal in the vicinity of the junction and the density of Cooper pair is 
exponentially descending together with the distance from the superconductor 
(see Fig. 2.2). The leakage of Cooper pairs into the metal can be understood 
as classical particle diffusion, where the pairing (coherence of Cooper pair) is 
gradually destroyed by scattering on other particles or lattice ions of metal. 
In the Ginzburg-Landau theory [89], the Cooper pair function is described by 


following macroscopic wave-function 


d(r) = V(r), (2.42) 


where N(r) is density of Cooper pairs and ¢(r) is their phase. 


Iv 


X 


Figure 2.2: The visualization of the Andreev/normal reflection at the metal- 
superconductor interface and the superconducting proximity effect (marked by 


red dashed line) for T < Tę, where £s denotes the coherence length. 
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The electron diffusion differs when particle transport is from metal to 
superconductor. In metal, a moving electron characterised by an energy Eko 
with a Bloch wave-vector k and spin o can cross the metal-superconductor 
boundary. When the energy of electron is smaller than the superconducting 
gap |e, — Ep| < A, where Ep denotes the Fermi energy of the metal, the 
electron cannot enter into a superconductor as a single particle (the normal 
reflection). In other possible process, an incoming electron can be bound 
to another electron with opposite momentum and spin, thus forming the 
Cooper pair which next propagates in the superconductor. This process, 
called the Andreev reflection, happens when the incoming electron with 
energy smaller than the energy gap is reflected as a hole with the opposite 
wave-vector —k and spin -o. The Andreev reflection involves the electron-hole 
transition process [90,91], the reflected hole moves along the same path with 
reverse direction as the incoming electron. The incoming electron crosses 
the metal-superconductor interface and combines into the Cooper pair with 
the remaining electron having the spin and momentum opposite to its own. 
In case of |e, — Ep| < A the reflection is dominated by Andeev process with 
high probability R4 ~ 1 and the normal reflection may be ignored. Otherwise, 
for energies larger than the gap, the electrons be normally reflected (see Fig. 
2.3a). 

The Andreev scattering can be described by the microscopic approach 


with the Bogoliubov-de Gennes equations [92,93]: 


Byla) = (gm u) ule) + A@)o(@), 
Ey(2) = ( so n) IEA. (2.43) 


2m dz? 


where u denotes the chemical potential, A(x) is a spatial pairing potential 


and y(x) is the 2-component wave-function 


Wart = (s (2.44) 
U 
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ee ea (2.45) 
defining the position in the space of the quasiparticles. By k* and k~ we denote 
the wave-vectors for electrons and holes, respectively. In such model, we can 
expect different Fermi momenta in the two materials, some elastic scattering 
in the electrodes and at the metal-superconductor interface. Assuming that 
there is no dissipation in the leads, the scattering at the phase boundary can 
be modelled by the delta function V (x)= Hê(x), with a potential barrier H. 
At the metal-superconductor interface the following boundary conditions can 
be assumed: the continuity of v at x =0, i.e. wy(0) =vs(0) and its spatial 
derivatives h/2m(%5(0) -v4(0)) = Hy(0). Using these assumptions, the few 
types of plane waves can be distinguished: Vinc an incident wave, a reflected 


wave refi and a transmitted wave Wransm, having two possible solutions [93]: 


Pinc = (es (2.46) 

ma 1 zk 
Wrefl = a( Jet + uf Je e, (2.47) 
WPiransm = de + a en (2.48) 


with 


Ak* = V2m[u + V E? - A2]"?, (2.49) 


hq* = V2m[p+ E]", (2.50) 
and the appropriate coefficients: a = uv/y, b = -(u? — v?)[(H/hupr)? + 
iHfhwp]/y, c = u(1-iH/fhwvr)fy, d = ivH/yhur, where y = u? + (u? - 
v?)(H/hup)?. Leaving aside the elastic scattering at the boundary of those 
two materials, the coefficients b and d are zero, which makes the Andreev 
reflection out of dissipation of quasiparticles. The amplitude coefficients u, v 


are energy dependent and they are defined in (2.29). 
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The Andreev scattering can appear in different heterostructures, e.g.: 
two N-S interfaces combined into a S-N-S heterostructure which form a 
Josephson junction [91], schematically presented in Fig. 2.3a. In a Josephson 
junction charge current flows without an applied voltage, this is known as a 
supercurrent. The intensity of the supercurrent is exponentially dependent 
on the correlation length of Cooper pairs, thickness of the normal metal layer 
and what is the most important the macroscopic phase difference between the 
superconductors © = ¢g; — s2. For the tunnelling processes the transmitted 


wave-function is phase dependent 
uer?si ee yei?s2 ee 
Wtransm = al Je” AE af je" (2.51) 
v u 


Due to the finite length of metallic material, the discrete levels are formed 
in the energy spectrum [92-95]. Satisfying the continuity conditions for 
wave-functions and their derivates at the material boundaries (at x = +L/2, 
where L is width of the normal metal) the dispersion relation for particle 


with energies below the superconducting gap |E| < A takes form 


c2ia(E) gi(k*-F")Lężić — 1 (252) 


where the energy-dependent phase factor is defined by a( E) = arccos(E/A). 
In the systems composed of high-density superconductors (with negligible 
spin-orbit interaction and without magnetic field), neglecting the elastic 
scattering and considering the particle energies much lesser than the Fermi 
energy |E| « Ep the eigenvalues can be calculated by using the relation 
h(k* -k')=2Efvp. It yields 


h 1 
Et = SE [2x (n + 5) + 3] (2.53) 
where n = 0,1,2,..., unless the levels approach A. In the case of Andreev 


limit |E| « Ep [96], the Andreev reflection probability reaches from R4 ~ 
T2/(2-T,)? at E = 0 upto Ry = 1 at E = A [97], where T, denotes the 


normal transmission coefficient. In point contact superconducting junctions, 
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Figure 2.3: Scheme of: a) Bound states in a short S-N-S junction with phase- 
dependent energy E(6) below the energy gap A, as result of Andreev reflection. b) 
Andreev bound states formed in S-QD-N junction, having essential role in transport 


spectroscopy [95]. 


where the superconductors are weakly connected, the contact region can be 
treated as the normal metal. Then, the metal region is much smaller than 
the coherence length of Cooper pairs Ly « és. For energies smaller that the 


energy gap |E| < A the relation in (2.52) has following solution 
a(E) = arccos(E/A) = +(6/2) (2.54) 
which implies a simple relation 
E =+Acos(¢/2). (2.55) 


Taking into account the elastic scattering, the discrete energy levels (bound 


states) are given by 
En = A[1- T,sin?(6/2)]"*, (2.56) 


with n = 1,2,...m the number of allowed modes and 7, denoting the normal 
transmission coefficient for each independent scattering-matrix eigenmode in 
the normal phase [98,99]. The scattering-matrix T, takes values from range 


(0,1). 
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In the metal-superconductor heterostructure, where the quantum dot is a 
barrier between the materials, the Andreev states can also appear (see Fig. 
2.3b). Emergence of the Andreev states in quantum dot strictly originates 
from the neighbouring superconductor. Their nature depends on the size 
of the pairing gap A and the Coulomb repulsion U between electrons [7]. 
For relatively large value of A, the multiple Andreev reflections give rise to 
Andreev bound states [100]. In the other limit, when the Coulomb repulsion U 
is dominant value in the system, the Andreev reflections are suppressed [101]. 
The influence of A and U on Andreev states will be discussed in details in 
Sec. 3.1. Emergence of the Andreev states will be systematically investigated 
in Sec. 3.4. Let us emphasize that such bound states are well detectable 
experimentally. Fig. 2.4a presents the recent results obtained for the Andreev 
in-gap states observed in differential conductance of current versus applied 
source-drain voltage and Fig. 2.4b shows their behaviour at the different 
temperatures. 

Moreover, the proximity effect can be found in the different type of 
systems, e.g. in superconductor-ferromagnet heterostructures [102], what 
gives opportunity to study the dynamical properties of Josephson nanojunction 
with ferromagnet impurities [103, 104]. 

In multiterminal devices based on quantum dots, the other types of 
Andreev’s scattering take place. For example, in a Cooper pair splitter device 
experimentally realised by Schónenberger's group [105] and theoretically 
investigated by [106-108] (schematically presented in Fig. 2.5), we can 
distinguish few types of scattering. The first possibility is ballistic transport 
(not presented in figure), which involves the flow of both electrons through the 
induced Andreev states on quantum dot, between the left and right electrodes. 
A second possibility for charge transport is the Direct Andreev Refflection 
(DAR), which occurs when an electron from the left electrode flows through 
the quantum dot on which it combines a pair with an electron from the second 


electrode (schematically presented in Fig. 2.5a). Then the paired electrons 
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flow into the superconductor. The third transport possibility is the Crossed 
Andreev Reflection (CAR), which involves the flow of an electron from one 
electrode to the quantum dot, on which the electron pairs up with an electron 
originating from the electron-hole transition process. In such case, the hole 
is reflected to the opposite electrode, what is schematically presented in Fig. 


2.5b. 
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Figure 2.4: The experimental evidence for Andreev bound states (ABS) in graphene 
quantum dot: a) the differential tunnelling conductance versus bias voltage, b) 
in-gap states observed for different temperatures: 0.26, 0.45, 0.67, 0.86, 1.25 and 
1.54 K [109]. 
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Figure 2.5: Schematic visualization of the multiterminal junction showing: a) Direct 


Andreev Refflection (DAR), b) Crossed Andreev Reflection (CAR) processes. 


We presented a few examples of devices based on quantum dots and 


CHAPTER 2. INTRODUCTION TO SUPERCONDUCTIVITY 34 


possible scenarios of charge transport. In this thesis, we consider only two- 
terminal nanojunctions with the single and double quantum dot. Their 


theoretical investigation will be presented in the next chapters. 


Chapter 3 


Dynamics of Andreev states in 


single quantum dot system 


This chapter describes dynamical properties of a single quantum dot placed 
between the metallic and the superconducting leads (N-QD-S), schematically 
presented in Fig. 3.1. We would like to investigate the time-dependent 
processes accompanying the formation of Andreev bound states in the system 
caused by the connection of circuit elements to each other. We identify 
their features and behavior depending on the system properties and initial 
conditions and determine the time necessary to reach equilibrium or steady- 
state by the system. Further, we investigate how changes in system settings 
affect the Andreev bound states. We analyse the evolution of the system after 
these changes, which reveals a number of new physical phenomena. In this 
chapter, we present results of time-dependent charge current and differential 
conductance, which could be experimentally examined. In what follows, we 
investigate the dynamics imposed by: transient effects (Sec. 3.4), rapid 
application of source-drain voltage (Sec. 3.5), quench in the orbital level (Sec. 
3.6) and quench in coupling to superconducting lead (Sec. 3.7). Moreover, 
this chapter includes the theoretical description of the mathematical method 


based on the Heisenberg equations of motion technique (Sec. 3.3), used to 


35 
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determine the time-evolution of the system. 


3.1 Model Hamiltonian 


The quantum dot (QD) coupled on one side to the metallic (N) and on the 
other side to the superconducting (S) reservoirs (Fig. 3.1) can be described 


by the second quantization Hamiltonian 
H = Hop + Hy + Hy-op+ Hs + Hs-gp. (3.1) 
The quantum dot part is given by 


Hop= > eadid, + U fyń, (3.2) 


where e, is its discrete energy level, dy (d,) is creation (annihilation) operator 
of the quantum dot electron with spin o = {t,{}. The second part of (3.2) is 
an interaction term describing the Coulomb repulsion U between the electrons 
with the opposite spin, where n, = did, is a particle number operator. 


The normal metallic lead is described as the free fermionic gas 
Hin > ZE (3.3) 
ko 


where summation is over the energy Eko of electrons with wave-vector k 
and spin a. In the system, the hybridization between the quantum dot and 


external lead is expressed by 
Hy_gp = 3 Vy (éld T h.c.) , (3.4) 
ko 


where Vą denotes the coupling of itinerant electrons with quantum dot. 


The superconducting effects are considered by the BCS-type Hamiltonian 
Hs= ) eu = >, (Acc) the), (3.5) 
qo q 


where Asc denotes the pairing energy gap. The hybridization between the 


quantum dot and superconducting lead is described by 


Hs-qD 5 > Va (cd, t h.c.) (3.6) 
qo 
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Figure 3.1: Illustration of single quantum dot placed between the normal (N) and 
the superconducting (S) reservoirs. The energy level of QD can be controlled by the 
additional gate Vg, U is the Coulomb repulsion potential and T'y/s are couplings 


to the external reservoirs [110]. 


with an appropriate coupling constant V}. Assuming this gap Asc as the 
largest energy scale in the system, we will focus on the processes in the 
subgap regime. Under such conditions the fermionic degrees of freedom can 


be integrated out, resulting in the low-energy physics described by |111, 112] 
A A A NE T X. os 
Hap + Hs F Hs-QD R Y ead} do +U Arh, t = (did, F h.c.) 5 (3.7) 


where Ts/2 plays a role of the effective pairing potential induced in quantum 
dot [61]. For simplicity, we consider the energy-independent coupling T's = 


2m X |Val° (€ - £4) in the wide-band limit approximation. 


3.2 Selected experimental realizations 


Presently available techniques for fabrication of nanostructures allow 
to create various devices based on the quantum dots. As an example, we 
present a few of such systems. The popular way to construct a nanojunction 
based on a single quantum dot is STM (scanning tunneling microscope). 
Using STM tip one can deposit either a single atom or a cluster of atoms 
on a surface of arbitrary material. This approach allows to control the 


coupling strength between the quantum dot and the surface by defining the 
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distance between them. It gives the possibility to study the bound states 
inside the superconducting energy gap by conductance spectroscopy [113]. 
Another example, Fig. 3.2a (discussed in Ref. [114]) presents multiterminal 
heterostructure composed of the carbon nanotube (CNT) partially covered 
by superconducting layer (S) of Pb. The carbon nanotube is connected 
on both sides to electrodes made of Pd (N;, Nə). In the system, SG, and 
SG are palladium controlling gates. The manufactured system presented in 
Fig. 3.2a can work as the simple normal metal-quantum dot-superconductor 
heterostructure and as the Cooper pair splitter [105], depending on the 
controlling gate potentials. Another example of N-QD-S system is presented 


a) b) 


AE: z eV 
'2 SG, 200 NM = 2 2 Back gate 


InAs/InP 
NW 


Figure 3.2: The experimental realisations of single quantum dot system. Panel a) 
presents the single quantum dot system based on the carbon nanotube |114]. Panel 
b) is schematic visualization of the single quantum dot system made of InAs/InP 


nanowire presented in panel c) [115]. 


in Fig. 3.2b,c [115]. The middle panel is a schematic visualization of the 
manufactured system presented in the right side panel. In this case, the system 
is made of the InAs/InP nanowire covered by the gold (N) and vanadium (S) 
leads. Under the nanowire the plunger-gate (pg) and superconductor-barrier 


gate (sg) are located. 
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3.3 Equations of motion approach 


For investigating the evolution of any observable O(t) we use the following 
equation of motion 
(OO), DOG Ty EOB: (3.8) 
dt h ot 
In particular, the time-dependent charge current flowing between QD and 
metallic lead jy,(t) can be found from the Heisenberg equation of motion 
for the particle number operator fiz, (t) = ĉl (t)êko(t). The current flowing 
in the superconducting reservoir js,(t) and occupancy of the quantum dot 
n,(t) = (a,(t)), can be indirectly calculated from the charge conservation law 
“no(t) = jno(t) + jso(t). 
We start from the momentum and spin dependent current jyko(t) defined 


as the time derivative of the statistical average for the particle number operator 


(za(t)) 


jnis) = =e (no(t) = -F (CH, fao (D). (3.9) 
Substituting the commutator 

[Fl fko(t)] = -Vći (t)dz(t) + Vel (t)êko (t) (3.10) 
into (3.9) we obtain 

jwio(t) = EW (Hd(t) - Vid (t) eet) (3.11) 


Subtracting the conjugated numbers z - z* = 2iIm(z) gives us the current 


formula for charge current 
; 2e + R 
JNko(t) = h Im (Vid! (t)Cxo(t))) é (3.12) 


The total charge current contributed by all electrons is obtained by summation 


over momenta k 


2 


ivel#) = Flm (FvO) (3.13) 
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At this step, the current formula depends on the quantum statistical 
averaging of operators (d;(t)ćko(t)), which is not trivial problem to solve. 
One possible route to obtain this relies on expanding the operator ĉko(t) into 
sum of the initial ć,,(0) and additional terms, emerging during evolution 
of the system. The equation of motion for annihilation operator ko (ż) = 
i[A,ć,,(t)] (for h = 1) yields 


Ów) + iEkoCko(t) = -i Vdo (t). (3.14) 


We can express (3.14) as a differential equation of the form 


dy(x) 


ae + f(x)x =y(x), (3.15) 


which has a general solution 
y(t) = y(O)e So Fed 5 fare! fedri (3.16) 
Thus, the explicit formula for time-dependent operator ć;,(t) takes the form 
Eo (t) = Ero (0) fo roA — i A “dt Vgert ieee"), (317) 


Substituting (3.17) into (3.13), we obtain the total current expressed by 


following equation 


volt) == D Im (WOBO (0) -i f Vee Vie Od(¢)), 

' (3.18) 
where V(t) = Vei fo Ero) denotes the time-dependent coupling to the 
metallic lead and V* (t')Vpo(t) = [Vp] 2e Jer <ro(t')dt" is a product of the time- 
dependent couplings at different times. In the wide-band limit approximation, 
assuming constant energy levels in time Eko (t) = Eko, the hybridization function 


V; (t')Vio(t) appearing in (3.18) simplifies to [116] 
- pt WY) 40 ert " Tr 
5 [Vee ret Jdt" _ Y [Ve] e e exo z => St), (3.19) 
k k 


where Ty = 27 P |V|? (€ — eg) is the coupling to metallic lead. 
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Finally, one obtains the charge current expressed in the general form 


jNno(t) = 2Im b Ve “rot (dt (t)Exo(0)) - = e(t) (3.20) 


To calculate the current jyo(t) we have to solve numerically a set of cou- 
pled differential equations for the expectation values: (A„(t)), (d;(t)d;(t)), 
(di (t)éx4(0)), (di (t)exy(0)), (dy (t)exy(O)) and (d,(t)é,)(0)). Explicit form of 
the individual functions is presented in Appendix A.1. Such method allows 
to investigate the time-dependent properties of the order parameter func- 
tion y(t) = (d,(t)d;(t)), which characterizes evolution of the on-dot electron 
pairing. 

For numerical purposes we use the Runge-Kutta of 4th order, assuming 


various initial conditions for occupancy of the quantum dot n,(t = 0). 


3.4 Transient effects and steady-state properties 


To demonstrate dynamics of the Andreev bound states we first consider 
their buildup after suddenly connecting the quantum dot to external reservoirs. 
Processes accompanying the evolution from the initial conditions to the 
equilibrium state, when the steady properties are established, can be regarded 
as transient effects. For computational reasons it is convenient to assume 
that initially, at t = 0, the quantum dot is disconnected from both external 
reservoirs and it is connected to the system constituents at t = 0+. Fig. 3.3 
presents numerical results obtained for currents jy+(t), js(t), occupancy of 
quantum dot n;(t) and the real part of the order parameter x(t). Imaginary 
par of x(t) is negligibly small in comparison to |Rey(t)|. We observe the 
evolution of the system (in absence of source-drain voltage) towards the 
steady state dependent on various initial quantum dot occupancy. The initial 
occupancy affects the transient currents and the charge accumulated in the 
quantum dot over time, as can be seen in their plots. If initially, the quantum 


dot is empty, we observe quantum oscillations due to the flow of Cooper 
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pairs back and forth between the quantum dot and the superconductor. In 
another case, when the quantum dot is initially occupied by a single electron, 
we observe gradual filling of quantum dot by electrons from the normal 
lead, because the Cooper pairs cannot flow onto the quantum dot. The 
relaxation processes originate from the coupling of the quantum dot to a 
continuum of electrons of the metallic lead. For the transient currents, the 
envelope functions is described by e""Nt [117]. We also observe that the initial 


occupancy of the quantum dot does not affect creation of an on-dot pairing. 
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Figure 3.3: The transient evolution of currents jn+(t), jst(t), occupancy m(t) and 
the order parameter Re x(t), imposed by a sudden coupling of the quantum dot 
to both external leads at t = 0* for various initial occupancy (n4(0),n,(0)). The 
results are obtained for zero source-drain voltage and system parameters: £q = 0, 


U =0, Ty =0.2, Tg =1. 


The oscillating character of evolution is reminiscent of the Rabi oscillations 


of two-level quantum system. At time t = 0, the quantum dot is characterized 
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by the single energy peak at ca. Then for the time t > 0, the quantum 
dot energy level evolves into a pair of Andreev bound states at energies 
E4= +/e2 + (T's/2)? . Those quantum oscillations are caused by fluctuations 
in the occupancy of these two-level system. In the transient evolution, the 
activation of quantum oscillations strictly depends on the initial conditions 
(for analytical reasoning see Ref. [117]). They are not observed if the initial 
configuration [n4(0),n,(0)] is singly occupied (1,0) or (0,1), as shown in Fig. 
3.3. 

Empirical observation of the Andreev states formation can be done by the 
differential conductance Gy(V,t) = z5jn;(t) of the charge current flowing 
between the quantum dot and the metallic lead, when the chemical potential 
of the normal lead uy is detuned from us by an applied source-drain voltage 
eV = uyn — us. We assume the superconductor to be grounded, jug =0. The 


voltage is imposed at t = 0. Fig. 3.4 shows the differential conductance Gy (in 
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Figure 3.4: Evolution of the differential conductance Gy (in units of 2e*/h) with 
respect to source-drain voltage V and time t, after a sudden coupling of the quantum 
dot to both external leads at t = 0*. For computation we used: eg=0,U=0, 
Ty = 0.05, Peat, 
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units of 2e?/h) with respect to the source-drain voltage V and time t. In this 
particular uncorrelated case of quantum dot system, where e, = 0, we observe 
formation of the Andreev bound states localized symmetrically with respect 
to the chemical potential, at energies E4 = +l s/2. Emergence of the Andreev 
peaks has an oscillating character with period 2r/E4, damped in time. Their 
damping is governed by the envelope function e""N* and originates form the 
coupling of the quantum dot to the continuum of free electrons of the metallic 
lead [117]. 

To calculate the conductance in the correlated quantum dot system by 
Heisenberg’s equations of motion technique we adapt mean-field approxi- 
mation, suitable for a weakly correlated system. Using the Hartree-Fock- 
Bogoliubov (HFB) approximation the Coulomb repulsion term in (3.2) can 


be written as 


didydid, = m4(t)dld, + ny(t)did, - n (t)r (t) 
+ x(tbdidj + x* (tdid; - x)’. (3.21) 


The Hartree-Fock term can be included in the renormalized quantum dot 
energy level €g = €q + Un_,(t) and the anomalous contribution rescales the 
effective pairing potential Is = Ts- 2U x(t) of the induced order parameter 
potential x(t) = (d;(t)d,(t)). 

Fig. 3.5 presents the results obtained for weakly correlated heterostructure. 
Here, the top and bottom panels show the results for Coulomb repulsion 
U = 0.25 and U =1, respectively. Both panels refer to particle-hole symmetry 
(ea = -U/2), when the quantum dot is half-filled in the steady state limit. 
Like in the uncorrelated case, the differential conductance Gy evolves towards 
its steady-state limit, exhibiting emergence of the bound states around E4 = 


+y/(€a+U/2) +(Ts/2)2. Relaxation processes are governed by function 


eT vt and the quantum oscillations have the period 27/F 4. 
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Figure 3.5: Evolution of the differential conductance Gy (in units of 2e?/h) with 
respect to voltage V and time t, after a sudden coupling of the quantum dot to both 
external leads at t > O treating the electron correlation effects by the mean-field 
approximation. In panel a) U = 0.25 and b) U = 1. The other model parameters 
are: Ty = 0.05, ITs = 1, eg = -U/2 [110]. 


3.5 Influence of suddenly varied source-drain voltage 


Transient evolution of the uncorrelated system has been characterized by 
the time-dependent observables shown in Fig. 3.3 and conductance shown 
in Fig. 3.4. To have a deeper insight into the time-dependent observables 
after the perturbation, let us discuss separately the plots of jya(t), Jsa(t) 
presented in Fig. 3.6. Until t = 30/Ts, we display the transient evolution 
caused by abrupt coupling of the quantum dot to the external leads of unbiased 
junction. The damped quantum oscillations are controlled by the coupling 


Ty to the normal lead [117,118], what is visible in the analytic formulas for 
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the observables: 
- Is -Tyt 
jso(t) = > sin (Tst) ^t, (3.22) 
T 
jno(t) = a sin (Tst + po) e TN*, (3.23) 


where gy denotes the phase shift between the currents. When the system 
reaches its equilibrium we apply source-drain voltage to the system at t = 
30/T s. We assume step-like change of the chemical potential Ly > uy + Aun 
at t = 30/T s. In Fig. 3.6, the evolution exhibits quantum beats appearing in 
the current jyo(t). 


0.2 i 


Aun = 10.0 
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Figure 3.6: The transient evolution of currents (in units of e/h), after a sudden 
coupling of the quantum dot to both external leads at t = 0* and their evolution 
after the jump of chemical potential level uy > uy + Aun at t > 30. The other 
model parameters: Ly = 0, Ty=0.2, Ts = 1, c4 = 0, U = 0 [119]. 


At this point, we would like to recall the essential properties of the 
quantum dot system without the presence of superconductor. It has been 
shown that the abrupt change of the source-drain voltage has significant 


influence on the current flowing through the quantum dot placed between 
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normal leads [120-122]. After the sudden application of the source-drain 
voltage, the time-dependent current structure shows the coherent oscillations 
with the frequency w = |UN -€al. 

In our system, the proximitized quantum dot is characterized by a pair 
of in-gap bound states with energies +T's/2, therefore the current jy,(t) 
is a superposition of two contributing currents with the frequencies wyją = 
jin +T 5/2, respectively. In consequence, the current jyo(t) is combination of 
the high frequency component w = uy and the low frequency beats with the 
time period T = 2r/T 5. Moreover, we notice that the current j5,(t) oscillates 
with only lower frequency. Difference between these two currents is caused by 
the time-dependent occupancy of the quantum dot, obeying the local charge 
conservation no(t) = jno(t) + jso(t). 

Analysis of the experimentally measured time-dependent current can 
indirectly give information about the coupling strength between the quantum 
dot and the superconducting reservoir, by determining the period of current 


oscillations T's = 2r/T. 


3.6 Quench in orbital level position 


Let us consider influence of the quench imposed on the quantum dot 
energy level. In this scenario, we abruptly change the quantum dot energy 
level eg > ea + Aeg at time t = 30/Ts (like previously, after the transient 
processes are completed) and observe the system’s response with completely 
different behaviour of evolution for all observables (Fig. 3.7). The post-quench 
evolution exhibits the oscillatory structure of all observables. The oscillations 
depend on the energy level of the quantum dot and the coupling to the 
superconducting lead T = 27/\/T'% + 4(Aeq)? [117]. In absence of the source- 
drain voltage (uyn = us), the oscillations gradually disappear over sufficiently 
long time and the currents jso, JNa vanish. In evolution after the sudden 


shift of quantum dot energy level, the envelope function has exponential form 
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Figure 3.7: The transient evolution of currents (in units of e/h), imposing a sudden 
coupling of the quantum dot to both external leads at t = 0* and their the post- 
quench evolution after the jump of the quantum dot energy level eg > €q + Aeg at 


t > 30. The other model parameters: uy = 0, Fy = 0.2, Ts = 1, eg = 0, U = 0 [119]. 


~ el xt (where t= t- 30/T s denotes the time measured from the quench of 


Ed). 


Influence of electron correlations 


We now consider the Coulomb repulsion U between two electrons and its 
influence on the post-quench evolution induced by the shift of the quantum 
dot energy level. In the system with the correlated quantum dot, the orbital 
energy level eg = -U /2 corresponds to the particle-hole symmetry point and 
it is convenient to take this value as the initial condition. Fig. 3.8 shows 
the evolution of physical observables after the quantum dot energy level is 


changed from e, = —U/2 to its final value £q indicated on the y-axis. 
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Figure 3.8: The time-dependent quantum dot's occupation n(t), current js(t) and 
the real part of (d,(t)d;(t)), after a sudden quench in quantum dot energy level 
from eg(t < 0) = -U/2 to eg(t > 0) = eg. The other model parameters: Ty = 0.05, 
Ts =4 and U = 1 [110]. 


In this case of strongly coupled quantum dot to the superconducting 
lead (Fs/U = 4), the evolution presents the damped oscillating behaviour 
of observables, similarly like without correlations. The amplitude of the 
oscillation increases with a greater difference between the quantum energy 
level before and after the quench |eg(t < 0) — ea(t > 0)|. The more the energy 
levels are shifted away from the levels before the change, the stronger we 
disturb the system, making its response more intense. This is better noticeable 
as we move away from the particle-hole symmetry point (indicated by the 
dashed line). When the occupancy of the quantum dot is far from half-filling, 
the Cooper pair can flow back and forth between the quantum dot and the 


superconducting lead, what is visible by the well pronounced oscillations in 
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the observables. Since QD spectrum consists of pairs of the Andreev bound 
states at +4 these dynamical processes are reminiscent of the two-level Rabi 
oscillations. As in the uncorrelated case, the damping process occurs with a 
relaxation rate T ~ 1/['\y. What is relevant, the perturbation in the quantum 
dot energy level affects the post-quench state of the system in the long time 
limit. For example, the occupancies tend to stationary value n(t > co) x 0.57, 
for the final QD energy ¢4/U = 0.5 and n(t > œ) » 1.23, for eq/U = -1, 
respectively. 

Another feature appearing in post-quench evolution is m-shift of the 
oscillations upon crossing the half-filling level eg = -U/2, marked by the 
dashed lines in Fig. 3.8. This phenomenon resembles the 0-7 transition 
reported for Josephson junctions under stationary conditions [123,124]. The 
oscillatory behaviour is also visible in the real part of the time-dependent 
order parameter y(t), which can be regarded as a quantitative measure of 
the on-dot pairing. In a short time after the quench, the value of Rey(t) 
oscillates with significant amplitude. 

In the case, when the quantum dot is weakly coupled to the supercon- 
ducting lead, 5/U = 1, we observe significant modifications of the oscillatory 
time-dependent quantities (see Fig. 3.9). In particular, the weaker the 
coupling strength I's enlarges the period of the oscillations of all quanti- 
ties. We observe that the final value of the quantum dot energy level after 
the quench has an extensive effect on the frequency. For a larger value of 
lea(t < 0)-eu(t > 0)| the shift of frequency to higher values is visible. Moreover, 
it is important to notice that the change of I's has no impact on amplitude 
range of the superconducting current js(t). Smaller value of I's is equivalent 
to the lesser value of the on-dot pairing potential. The change of I's affects 
the final values of the other quantities. For example, n(t > oo) = 0.2 for the 
quench to e4/U = 0.5 and n(t > oo) x 1.55 for eg/U = —1, exceed the range of 
the occupation number in Fig. 3.8. 


In addition, to complete our investigation of the system evolution after 
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Figure 3.9: The time-dependent quantum dot's occupation n(t), current js(t) and 
the real part of (d,(t)d;(t)), after a quench in quantum dot's energy level from 
ealt < 0) =-U/2 to eg(t > 0) = eg. The other model parameters: Ty = 0.05, I's = 1, 
U = 1 [110]. 


shifting the quantum dot energy level we present the differential conductance 
Gy diagram in Fig. 3.10. In this plot, we demonstrate the evolution of 
Andreev states and their steady state positions after the shift of the energy 
level eg = U/2 > -U /2. We visualize the post-quench evolution in the wide 
spectrum of the source-drain voltage. After the sudden jump of the quantum 
dot energy level eg = U/2 > -U/2 we observe a similar response of the 
system like in previous cases. In particular, the quantum oscillations of the 
differential conductance are damped in time and the Andreev bound states 
emerge with the same spectral weight located with respect to the chemical 


potential uyn = eV. 
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Figure 3.10: The post-quench evolution of the differential conductance Gy (in 
units of 2e?/h) with respect to voltage V and time t, after a sudden change of the 
quantum dot's energy level eg = U/2 > -U/2 at t=0. The other model parameters: 
Ty = 0.05, Ts=2,U=1 [110]. 


3.7 Response to sudden coupling [ls 


For further understanding of the dynamics of the perturbed quantum dot 
system, we investigate another type of quench. In this section we look at 
evolution after a rapid coupling of the quantum dot to the superconducting 
lead. Let us remind the stationary solution for the system in the limit of 
Iy = 0 and A > œ, where the model parameters eg, U and I's determine 
character of QD ground state. The quantum dot ground state can exist in 
the singly occupied |ø) or in the BCS-type u|0) - v|tl) configuration [125]. 


The border line between the doublet and singlet state occurs at 
Acree = U?, (3.24) 


where 67 = €a + U/2. Crossing from one to another configuration is called 
a quantum phase transition |61] and plays an essential role in the interplay 


between the on-dot pairing and the correlation effects. Here, Fig. 3.11 shows 
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Figure 3.11: The time-dependent quantum dot's occupation n(t), current js(t) 
and the real part of (d,(t)d;(t)), after a quench in the coupling strength Ts(t) 
from zero to its final value I's. The other model parameters: eg = -U/2-U/20, 
Ty = 0.05 [110]. 


the evolution of physical quantities after abrupt coupling of the quantum dot 
to superconducting lead, 0 > T's, for the half-filled quantum dot (eg = -U /2) 
and weak coupling to the metallic reservoir, ly = 0.05. As before, we observe 
the oscillating behaviour of the physical quantities in time, which is the result 
of the Cooper pairs leaking into the quantum dot and creating the Andreev 
states. In the system, the quantum phase transition occurs at I's =U when 
the quantum dot level y= -U/2. In Fig. 3.11, by the white dashed line, we 
marked the theoretical border between the singly occupied and BCS-type 
ground state configurations. The quantum phase transition is clearly visible 


in change of Rey(t) in Fig. 3.11, where we observe development of the order 
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parameter after coupling quantum dot to the superconducting lead. In the 
doublet region (Ts < U) the order parameter has negligible values. In the 
same region for js(t) (middle panel), we do not observe significant charge 
flow due to the dominant Coulomb repulsion. The singlet region appears for 
Is > U, where the system relaxes to the BCS-type ground state through a 
sequence of damped quantum oscillations with significant amplitude right 
after the quench. For stronger coupling I's we observe faster oscillations. The 
boundary between doublet and singlet states is not sharp. We observe a 
crossover rather than a transition between these configurations caused by the 
coupling Ty # 0. 

Quench in the coupling to the superconducting lead can be analysed in 
the differential conductance plots. Fig. 3.12 presents the results obtained 
for I's =0 > 2. For this specific scenario, we consider the half-filled quantum 


dot eg = -U/2 and set the Coulomb potential U = 1. Before the quench, the 
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Figure 3.12: The post-quench evolution of the differential conductance Gy (in 
units of 2e?/h) with respect to voltage V and time t, imposed by a quench in 
hybridization strength T's = 0 > 2 at t=0. The other model parameters: eg = -U /2, 
Ty = 0.05, U = 1 [110]. 
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correlated quantum dot is characterized by two quasiparticle peaks at energies 
ej and €g+U. The rapid coupling to the superconducting lead induces the 
proximity effect which entails the creation of new states at energies +E 4. 
Additionally, the energies of new-created Andreev states, in the limit ofTy=0 
are Ea ~ (ea + U/2)? + (T's/2)2. In the limit of Ty =0 the border between 
the singlet and the doublet configuration is well specified by (3.24). From the 
cognitive point of view, we can ask what is the response of the system when the 
coupling to the superconducting lead drives the system from one to another 
configuration. In Fig. 3.13 we present two cases, the first when we observe the 
transition from the doublet to singlet configuration by increasing the coupling 
I's (upper panel) and the second situation when change of the value of I's 
induces a conversion in the opposite direction (lower panel). Starting from 
the initial conditions, at t > 0, the quantum dot is coupled to both reservoirs 
and its energy level is detuned to the half-filled occupation eg = -U/2, where 
the upper panel shows strong coupling to superconducting lead with I's > U 
and lower panel a weak coupling case I's < U. Until the time t = 5/2P'y, both 
panels show the transient effects and formation of the doublet /singlet state 
on the quantum dot. When the system reaches equilibrium, the quantum 
oscillations of Gy will no longer occur. We abruptly reverse these couplings 
Ts at t= 5/20 y. This shift of T's value caused transitions from the doublet- 
to-singlet (panel a) and from the singlet-to-doublet (panel b), respectively. 
Regardless of the type of transition, we observe the quantum oscillations with 
period T = 27/E4, where E4 ~ (ea + U/2)? + (Fs/2)2 are Andreev bound 


states in the limit of Fy = 0. Moreover, the post-quench evolution is not 


completely identical to the transients of Gy. Differences have origin in the 


slightly varied initial configurations of the system before the quench. 
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Figure 3.13: The evolution of the differential conductance Gy (in units of 2e*/h) 
with respect to voltage V and time t, obtained across the doublet-singlet transition 
due to the quench of Ts performed at t = W The upper panel presents transition 
from doublet to singlet: I's = 0.6 > 1.4. The lower panel present transition from 
singlet to doublet: Ts = 1.4 > 0.6. The other model parameters: eg = —U/2, 
Tx = 0.05, U=1 [110]. 


3.8 Summary 


In this chapter, we have analysed the transient effects originating from 
linking the system components together. We have shown that the activation 
of quantum oscillations in transient charge currents strongly depends on 
the initial conditions. They are not observed if the quantum dot is initially 
occupied by a single electron. Otherwise, when the quantum dot is initially 


empty, the evolution to the steady state of the system exhibits quantum 


oscillations with the period T = 2r/Ey, where E4 = +y/ (sq + U/2)° + (T 5/2)? 
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is energy of the Andreev bound states. 

We noticed, that in the evolution of the system (uncorrelated quantum 
dot) after the sudden change of source-drain voltage the beats in the current 
JNa(t) appear. The current jyo(t) is a superposition of two contributing 
currents with the frequencies wyje = uy + V's/2, respectively. 

For the uncorrelated quantum dot, the evolution after the shift of orbital 
level position exhibits the oscillatory structure of observables. The oscilla- 
tions depend on the energy level of the quantum dot and the coupling to 
the superconducting lead T = 27/\/T'2 + 4(Aeq)?. Investigating the system 
including the correlations, we have shown the r-shift in the oscillations upon 
crossing the half-filling level, which resembles the 0 — m transition found for 
Josephson junctions in stationary conditions. 

Finally, we also investigated the evolution of the system from its singlet 
configuration to doublet state and vice versa by changing the coupling strength 
I's. This dynamic change of configuration activates the quantum oscillations 
with the period T = 27/E 4 regardless of the direction of the transition, where 
Eą = +4/(€4+U/2) + (Ts/2)2. 

Independently of the initial state, the relaxation of the post-quench pro- 


cesses are governed by function e^t, 


Chapter 4 
Periodically driven quantum dot 


In this chapter, we consider another dynamical scenario for N-QD-S 
heterostructure, where the quantum dot energy level is periodically varied by 
external fields. We investigate the influence of external fields on the Andreev 
bound states, the charge current averaged over a period of oscillations and the 
differential conductance. The periodically driven quantum systems constitute 
a broad field of physics, and we would like to explore the phenomena that 
appear in our N-QD-S nanojunction under such conditions. The suitable 
method to study systems described by periodic Hamiltonian H(t) = H(t+T), 
with period T = 27/w, is Floquet theory [126,127]. We will present the basics 
of Floquet theory and its application to the considered quantum system. 
Results presented in this chapter were obtained by combining the Floquet 


theory with Green’s function method. 


4.1 Floquet theory 


Let us start by presenting essential ingredients of the Floquet approach. 
For the time periodic Hamiltonian the solution of time-dependent Schródinger 
equation 


hS) = AOWA) (41) 


59 
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takes the special form called Floquet state |w,,(t)), where 


lb (t)) = fun (t) jet” (4.2) 


is time-periodic mode |u,,(t)) = |u(t + T)) (called micromotion of the system) 
and e, is the corresponding quasienergy. These Floquet states (4.2) are 


eigenstates of the time-evolution operator over one driving period 


U (to + T, to) |Yn(to)) = ern (to)). (4.3) 


The eigenvalue e~*7/* (the phase factor) is independent of to, at which the 
evolution over one time period starts. The time-dependent Floquet states 
\w,(t)) can be computed by acting with the time-evolution operator on the 
Floquet state in specific time to: |d(t)) = U(t, to) |Yn(to)). The time-evolution 


operator can be written in the form 
U (tą, t1) EA las, (t2))(u(t1)|, (4.4) 


because for every Floquet state one can choose a complete orthonormal basis 


at arbitrary time t. Then the evolution of a state |w(t)) can be expressed as 
kb(t)) = Dene enn (t), (4.5) 


where c, denote the time-independent coefficients, c, = (Un(to)|W(to)). 

The dynamics of the quantum system over times longer than the single 
driving period can be studied in a stroboscopic manner. In this approach the 
evolution is studied in the steps of the driving period T. This stroboscopic 


time evolution is simply expressed by 
U (to + T, to) = e o” (4.6) 
and depends of the Hamiltonian HE represented in the Floquet space 


He = > Enlun(to)) (Un (to)|. (4.7) 
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It is important to notice that the quasienergies e, and the Floquet modes 
Un(t) are not uniquely defined. The shift of quasienergy e, by energy hw 
does not affect the phase factor e~*»7/". Due to this fact, the quasienergy 
En can be chosen arbitrarily and their shift by energy hw can be performed 


multiple times, then the quasienergy can be marked by additional index 
Enm = En + Mhw. (4.8) 


Next, the corresponding Floquet mode takes a form |uy,(t)) = e”™*|u,(t)). 


At this step, the solution of Schrédinger equation has a general form 
lWn(t)) = SSB mę) = SEM (tN (4.9) 


Substituting (4.9) into the time dependent Schródinger equation (4.1) we 


obtain 


LO 3 in) RARONSSADE JI (4.10) 


in an extended Hilbert space F = H © Lr. The Floquet space F is a product 
of the spatial space H for quantum systems and the space Lr of square- 
integrable functions that are periodic in time. In the extended space F, the 


scalar product of states can be averaged over a time period as 


Uuo) == f dttd), (4.11) 
where double braket notation is used for elements of F, corresponding to state 
at time t and its spatial coordinates in H. The left side of (4.10) is called 
a quasienergy operator Q(t) =H (t) - ih. To distinguish the quasienergy 
operator in the Floquet space F and Hilbert space H, let us denote the 
quasienergy operator by Q if it acts on F. Then, the quasienergy equation 
can be written as Q|unm)) = Enmltnm)) in F space. 

In the Floquet space F the set of orthonormal basis states |am)) can be 


constructed as a combination of the orthonormal basis states |a) of H and 


the set of time-periodic functions e”* as |am(t)) = |a)e'"e'. In this basis 
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the quasienergy operator Q takes block structure 
INA 1 > -im'wtj I Fj „d imwt 
((a'm'iQlam)) = — dte Ha A(t) - ih lade 
T Jo dt 
= (a'|Hm_nla) + ómmówamhw (4.12) 


where Hy, = Lej: dte~™t H(t) = Hi, is the Fourier transform of H(t). This 
block structure of Q is also dubbed Floquet matrix Hamiltonian HE „ which 


has the following form 


F 1 d inwt A —imwt 
Hez T af dte"*"'Q(t)e 
0 
= Hm — MAWS ym. (4.13) 


It acquires the block structure 


A A 


Ho + lhw A Hə 
HE 2 |... A, A H, M ie (4.14) 
Ho Ay Ho - lhw 


where the diagonal elements are shifted by hw which are related with absorp- 


tion or emission of "energy quanta" from the environment. 


4.2 Solution for time-periodic Hamiltonian 


The Green’s function formalism (Sec. A.3) and Floquet theory (Sec. 4.1) 
can be combined and used to describe the quantum system driven by periodic 
external fields. For example, the single quantum dot Hamiltonian presented 


in (3.1) can be simply converted to the time-dependent form 
Ê = Hop(t) + Hy + Hy-op + Hs + Hs-qp, (4.15) 


with an oscillating quantum dot energy level Hgp(t) = X„ea(t)dyd,, where 


Ealt) = £9 + Acos(wt) and £9 is a stationary quantum dot energy level. 
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For the proximitized quantum dot it is convenient to introduce the matrix 


Green’s function Ggp(t,t’) (Nambu representation) 


RZE ~i(To{d;(t); di(t))) - i(Totd,(t);d,(t'))) (4.16) 
QO EIA w ity Ka RENE ZMOD E 
-i(totdj(t);di(t)))  -i(Totd,(t);d,(1'))) 
whose elements come from the general definition for the contour ordered 


Green's function 
a g(t t) = -i(To{ A(t); B(t’)}). (4.17) 
From the equation of motion for Green’s function Ge at, © t’) we obtain 
iGS (tat) = Belt boo = i(Fo{iOda(t)s dy, E) = 
b(t -t') Soar + cat) Ge, at, (t,t') + 2 VG rai, (t,t) 
D a): (4.18) 


where the mixed Green’s functions 


Ge odt, (bt) = -ifo {êsa (t); dht), (4.19) 
Gi. at (t, t") = —i(TotENko, (t); dL (t’)}), (4.20) 


describe the particle propagation between the superconducting /normal elec- 
trode and the quantum dot, respectively. For uncorrelated quantum dot the 


equation (4.18) simplifies to 
(ið; — ea(t)) Gy, at, (t,t) = 6-(t -t’) dco" + 2 E a (t,t') 
+ >. VAG ernest, (t,t). (4.21) 
The contour-ordered Green’s function can be rewritten as 
Chg (Ot) = boog, (+ EV fi, ENE (TE) 


OMe J dro ONG, „z, (MT). (4.22) 


CHAPTER 4. PERIODICALLY DRIVEN QUANTUM DOT 64 


where 9, it (t,t’) is bare Green’s function of isolated quantum dot. Determi- 

nation of the mixed Green’s function is discussed in Appendix A.8. At this 

step, let us consider the particle propagator Ge 46t) To solve (4.22) we 
DĄ 


need the following functions 


zat DV f drat (DG (rt), (423) 
DĄ 


CNkblN 


c NN c c I 
Ge (it) = -Va J. drad a ETG) (Tt) 
c c I 
+ Jw, (t, T)VGGA, jt ). (4.24) 
Substituting (4.23) and (4.24) into (4.22), we obtain 
GP) =I gE) 
faf, dr’ I, al a ane AG T,T)Gź Aor) 
1 I „I 
+ far far TEA IG at (7 b) 
c c ! c ! 1 
š J dr JL drag (ŁZE, (TIC, (TL) (428) 
where the contour-ordered self-energies are defined as 


a T')= 2 KE EJ; (4.26) 


Yiv, ch Tenine 

A ch, (7,7) = 2 Va 9, a (TT) (4.27) 
spre S-qt? Sql 

> ch AG JE L Va I “ad „© T’). (4.28) 


In the same manner, we can s the other element of the Green's 


function (4.16) in the Nambu representation 
c I _ rc c ! c r „I 
Gey (bt) = ie dr p dr'g „(DZE „(TYG (TP) 
! „Cc ! I „I 
+ |dr | dig ETB, 4 (> TYŚ Ae a) 
! „Cc c ! c ! 4I 
+ J. dr J dr gą „(ZĘ a (TGS ait). (4.20) 
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The solution for diagonal elements of (4.16) is sufficient to find other related 
functions. According to the Langreth’s theorem (A.57), the contour-ordered 


casual functions can be expressed in the terms of retarded and advanced ones 


Góry KE) =o abt) 
+f dr | dr'g, p PDZ. aż T)Ga are) 
n 8 dT > tg, CTE a 10 TEA KŻ 
+ jA dT arg at (SO C T)G a Tot) (4.30) 


Geet = fdr far’ ENEE ret) 
+f dr f drige EEF g (TYG (TL) 
f dr | dr RESIZE (4.31) 


At this step, we can perform the Fourier transform of the Green's function 


oo 1 T 1 R 7 
Sits) s i, a= |, dteiletnw)t-ile+mw)t Ft), (4.32) 
—00 0 
where the contour integral of the functions C(t,t') = J„drA(t,r)B(r,t') 


is equal to Chm(€) = Xx Ank(€)Bkm(e) in Fourier space (the proof of this 
property is included in Appendix A.5). Finally, we obtain the Green’s 


functions expressed by the Fourier coefficients 
r,a 
Gy d; aN e) = a. CJĘ ae e) + 29, „dł; NIĆ AE e)G; A “dt mte) 
r,a r,a 
+ ZE até OAZA dt sim (£) 


* Iha E e)> A) i au JG diam (£), (4.33) 


Gad. nm) = daw 07, CNT mc IG alam (=) 
EDI. sd, MODE chuck, at OG dll (e) 
k,l 


a EE aa E (4.34) 
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At this step, we can recast (4.33) and (4.34) into the Dyson equation in 


matrix form 
GoD nm() = SQDnm(E) a 2. BOD,nk (E)Z rr (EVEGp mE): (4.35) 
k,l 


where Miim(e) = ER n (€) + UC" (€) is the total self-energy of the system. 


N,nm S,nm 


The individual self-energies (4.39) are defined as follows 


Ty 
See | mae (4.36) 
N,nm E 0 Ty | 
2 
Na. „(EJ=F AGA o AB â (4.37) 
ane (€+40+)2+AŻ]|-A., E| M 


where a(ć) = O(Asc - |El) + isgn(£)O(|<| - Asc) and £ = € + nw. More specific 
calculations are presented in Appendix A.9. 

According to Ref. [128], the Dyson equation (4.35) can be solved when the 
self-energies of the system are diagonal in Floquet space. Then the solution 


of the matrix equation (4.35) becomes 


r,a r,a ra -1 r,a 
[G$p(e)] = [I- g$p(e)2"*(e)] [egp(e)], (4.38) 
which can be explicitly shown as 


-1 


Y-1-1 Y-1,0 7-11 
[Gae] — 0-1 wo Wa - | [eaae] (4.39) 
Y1-1 710 ıı 


where Ynm = Onml — SOD nm?1mm and I denotes the identity matrix. 
The bare Green's function SOD nm of isolated quantum dot can be diago- 
nalized with respect to the Floquet coordinates n, m by appropriate unitary 


matrix A [129] 


2, Am(8GD)nmAmi = (Dob) - (4.40) 
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In this basis, the diagonal retarded and advanced Green’s functions are 


expressed by 
(Dop)u = (€+ lw £10" )I - £402, (4.41) 


where o, is the Pauli matrix. Making back projection, the isolated bare 


Green’s function can be written as follows 
SÓDnm(5) = © AuDobu(E)Am- (4.42) 
7 


Since we consider oscillating energy level eg(t), the matrix operator A, is 


expressed by the Bessel functions of a first kind (see Eqs. (A.82)-(A.95)) 


A 
Annis Inm (3) : | (4.43) 


0 gaa 


Then the matrix representation of the bare Green's function ((A.96)-(A.97)) 
simplifies to [66] 


Jn-( Ale) Imi Al) 0 
r,a e+i0t+lw-e 
BOD,nm (E) = > 0 . In_1(-A/w) Jm—1(-A/w) . (4.44) 
í e+i0*+lw+ey 


In the case of time-periodic driving of the system, the effective quantum 


dot spectrum is determined by 


(Pa(e)) = -+ Im [Gópo(e +10*)],,. (4.45) 


which is equivalent to averaging over a single period T [128]. In general, this 
spectral representation function must be positive [130]. 

To calculate time-dependent charge current we use general formula intro- 
duced in (3.11) for particle transport into the metallic lead [65,131,132]. The 


charge current equation can be written in the form 


Cho dl, 


Inolt) = p D Ve Re [He (0: GO] = TV Refa; (£,6)|, (4.46) 


CHAPTER 4. PERIODICALLY DRIVEN QUANTUM DOT 68 


where the lesser Green’s function can be expressed as a convolution of Green’s 


function for quantum dot and bare Green’s function 


CNko;Ć 


Gs (6) =V f dr c= CoE (itt GS a OE MM JG). 


(4.47) 


Substituting the lesser Green's function to (4.46) we obtain the Landauer’s 


formula 


M 2e T a 
nett) = p frarRelor „(DES 9 D+ p ODE, g, 08) 
(4.48) 


Using the Nambu representation, the Andreev current (flowing through the 


normal lead) can be expressed in the matrix form 


jn) = 5 f dr Re[Gaplt, EN, t) + Ghot DER Oia (449) 


where the diagonal elements 11 and 22 correspond to the particle and hole 
terms, respectively. Transforming (4.46) into Fourier space, it is possible to 


express the time-dependent current in Floquet space 


wie de —UnN—-p)w r 
imd D | g PA OP Go amle) ml) 
himo w/ 


+G Op nm (E) ZN mp (€) | 11-22: (4.50) 


Assuming that the indices n, p are equal, we obtain the average current over 


single period T 


P de p r a 
Cn) = FE fg 27 ROG AD le) Ekle) + Spr (EB, le) 


(4.51) 


Here, in (4.50) and (4.51) the lesser Green’s function of a single quantum dot 
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in Nambu-Floquet space is equal to 


[Góonm(E)], 2 Od nk lm [Ekla LG Qn tm jaw 


+ [Gop nk] uE iE) 121 GO tm lav 
+ [Gop nk] ul File) 211 Gop tm div 


a [Gop nel pal Ee (€) J221G Od 1m lav}, (4.52) 


where X< (e) = U5,(€)+ i$ (€) is sum of the lesser self-energies of the reservoirs. 
The lesser self-energy for single electrode is the result of the subtraction of 


advanced and retarded functions 


HY nm(E) = (= en oe - SN abe SK SE; e) (4.53) 
HS nm(e) = (SS an (E) - ZS am(E)) fs(e + nw), (4.54) 


multiplied by the Fermi-Dirac distribution function 


fNjs(z) = Uf [e@ erika? + 1], (4.55) 


4.3 Numerical results 


The method applied in this chapter, based on the Green’s function and 
Floquet approach, allows to calculate the time dependent observables and 
observables averaged over the single period. In this method the system is 
adiabatically switched on and the time-dependent functions describing the 
system become cyclic with the same period as the drive. All results in this 
chapter are obtained for the temperature limit T > 0. As an example of the 
time-dependent currents we present Fig. 4.1. In this figure, we show the 
charge current Iy(t) = jn;(t) for different values of oscillation amplitudes 
in regime of the superconducting atomic limit Ase > oo. For the system 


described by Hamiltonian (4.15) with applied source-drain voltage V = lw, we 
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observe oscillatory character of the plots in single period window of the drive. 
At this stage, we would like to highlight the complex nature of these charts. 
For a larger amplitude, a more wrinkled structure of the plot can be seen. The 
reason for this phenomenon is that there are more eigenenergies in a transport 
window. Localization and behaviour of those eigenenergies will be revealed 
in the further part of this chapter. For the specific time window, the time- 


dependent currents obey the symmetry relation Iy(-V,t) = -In(V,t+ 7/2). 


0.075 
0.050 
= 0.025 
RR 
= 0.000} 
Rs 
—0.025 — A= bia 
— A=23/w 
—0.050 — Aj 
0 T/2 T 


t 


Figure 4.1: The time-dependent currents for the periodically driven eg(t) obtained 
for: V = lw, Ty = 0.lw, ITs = lw, U = 0. The results were obtained from the 


numerical solution of (4.50). 


One of the convenient ways to study the periodically driven system is to 
compute the observables averaged over a period. Using the Landauer formula 
in (4.51), we computed the time average current (Jy) and differential conduc- 
tance (Gy) averaged over one driving period versus the source-drain voltage 
V applied across the junction, considering the in-gap states by assuming 
A, > œ. In Fig. 4.2, we can notice that in the absence of source-drain 
voltage the current (Ij) is zero, because the incoming and outgoing charge 


flows cancel each other. For larger source-drain voltage V, we observe an 
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increasing value and plateau regions of (Jy). Using standard definition of the 
differential conductance (Gy) = m we reveal the energy states appearing 
in the driven system. We notice that the conductance peaks perfectly coincide 
with the slope parts of charge current plots, which exhibit the positions of 
eigenenergies. The differential conductance (Gy) shows that the rapid growth 
of the averaged current occurs if the voltage V comes close to the eigenenergies 


of the system. 


eV/w 


Figure 4.2: The averaged current (Jy) and differential conductance (Gy) with 
respect to voltage V obtained for the periodically driven eg(t), using: Ty = 0.1w, 
I'g= le, U = 0, 


To gain a broader view of the quasiparticle eigenenergies, let us present 
results for the spectral function (pg(<)) defined by (4.45). In the time-periodic 
system described by Hamiltonian H(t) = H(t+T') the spectral function is 
equivalent to averaging over a period T. In Fig. 4.3 and Fig. 4.4 we present 


the spectral function for N-QD-N and N-QD-S heterostructure, respectively. 
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Figure 4.3: The averaged spectral function of QD versus the amplitude of driving 


field. Results are obtained using: Ty = 0.lw, Tg = 0, U = 0 assuming ey(t) = 


Acos(wt). 
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Figure 4.4: The same as in Fig. 4.3 but for the N-QD-S heterostructure, assuming 
Tg = lw. 
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The results for the normal junction will be a reference point. The case of 
N-QD-N heterojunction in Fig. 4.3 shows the spectral function dependent on 
energy e and driving amplitude A, which include a series of equidistant states 
at energies £9 + nw, where n is an integer number. Those states are called the 
photo-assisted tunnelling states (PAT), harmonic states or Floquet states for 
the driven system [65-67]. We observe the emergence and extinction of these 
states and their greater number with larger amplitude of the drive. The zero 
harmonic peak vanishes for ratio A/w = 2.405 which corresponds to the zero 
value of Bessel function Jo(A/w). We see the finite broadening of the peaks 
which comes from the coupling to the metallic reservoirs. 

For N-QD-S heterostructure, such spectral function dependent on energy 
level e and driving amplitude A is presented in Fig. 4.4. In contrast to the 
N-QD-N junction, in this particular case, we observe the splitting of the zero 
and higher-order harmonics. In the stationary case (A = 0), the quasiparticle 
peaks in the subgap spectrum coincide with two Andreev bound states at 
energies +/(e9)? + (T 5/2)? [61] broadened by Ty. For small amplitudes, 
the evolution of the quasiparticle energies shows that the branches of the 
zero-harmonic mode split. Simultaneously new higher-order harmonic states 
appear, shifted by energy nw. We observe that the zero-harmonic states 
linearly branch out and they are headed to the horizontal lines corresponding 
to energies nw. One pair of the zero-harmonic branch crosses each other 
around A = 1.lw on the line corresponding to the zero energy. The second pair 
of the zero-harmonic state evolves to the horizontal lines for n = 1 and n = -1, 
where they interlace with the higher-order harmonics. For larger amplitudes, 
we observe overlapping and interlacing of the states between neighboring 
harmonics. Together with the larger amplitude the distance between mixed 
branches and their intensity slightly decrease. 

Comparing the differential conductance (Gy) in Fig. 4.2 with the quasi- 
particle spectrum for corresponding amplitudes, we see that the additional 


states in charge transport perfectly coincide with the positions of the harmonic 


CHAPTER 4. PERIODICALLY DRIVEN QUANTUM DOT 74 


states. The slightly smaller splitting of mixed branches along n = 0 comes 
together with the lesser expectation value of the on-dot pairing potential 
(dydyr = - JB Z iaie (4.56) 
averaged over a period T. Its dependence on the amplitude A is presented in 
Fig. 4.5. We observe that this induced order parameter is correlated with 
the amount of spectral weight of the zeroth-order harmonic. The proximity 
induced order parameter vanishes for such amplitude where the zero-level 


harmonic states lose their spectral weights. 


Order parameter 


4 
A/w 
Figure 4.5: a) Expectation value of the proximity induced order parameter (d dą) 
versus the amplitude A, obtained for several values of Ig, as indicated. Using: 


ed = (, [y/w =0.1. b) Zoom of the plot shown in panel a. 


To investigate the influence of the driving frequency, in Fig. 4.6, we present 
the time-averaged density of states versus w/l's. The results are obtained for 
the amplitude A = 2.21's. This quasiparticle spectrum illustrates one pair of 
the Andreev quasiparticles (zero mode) with the higher-order harmonic states 
and their linear evolution along nw. The higher-order states are aligned along 
eV + nw lines and disappear for large w. In the high frequency limit w > o, 
the numerical results agree with the analytical solution of Magnus expansion 


for the proximitized system Hamiltonian (discussed in the Sec. 4.4). 


=J 
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Figure 4.6: The quantum dot's spectral function of N-QD-S heterostructure versus 


the driving frequency w. Results are obtained for A=2.21'5, Fy = 0.1I 5, U =0. 


Influence of the pairing gap Asc 


To investigate a more realistic situation, we would like to focus on a case 
when the superconducting energy gap A,, takes a finite value. Usually, the 
value of the energy gap is a few or fractions of meV. Fig. 4.7 presents the 
quasiparticle spectral function with respect to the amplitude A, obtained 
for Ase = 0.5w. In this plot, we observe that the quasiparticle function does 
not follow the main properties of the energy spectrum displayed in Fig. 4.4. 
We can notice that the higher-order harmonics (|n| # 0) are pushed out of 
the superconducting gap and their splitting is reduced in comparison to the 
results in the limit Ase > oo. This behaviour depends on the efficiency of 
the pairing which gradually decreases with increasing amplitude, due to the 
larger influence of the continuum of electron states above the superconductor 
energy gap. Raising the QD energy level above the superconducting energy 
gap causes a partial leakage of the spectral weight towards the in-gap regime. 
A feature of this phenomenon is a continuous background which appears in 


the spectral function, what is visible in Figs 4.7 and 4.8. This continuous 
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background corresponds to incoherent subgap states. 


3 
2.5 
2 
2.0 
1 
3 1.5 
> 0 
if 1.0 
=; 0.5 
=3 
0 1 2 3 4 5 


A/w 


Figure 4.7: The quantum dot’s spectral function of N-QD-S heterostructure with 
respect to the driving amplitude A. The results are obtained for Ase = 0.5w, 
Ty =0.1lw, 9 =1w, U =0. 
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Figure 4.8: The quantum dot's spectral function of N-QD-S heterostructure with 
respect to the driving frequency w, using: Asc = 0.5T' s, A=2.21 s, Ty=0.1[s,U=0. 
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To illustrate the influence of driving frequency in the case of the finite 
value of the superconducting energy gap, we present Fig. 4.8. This figure 
displays the distribution of the spectral weight between the multiple harmonics, 
presents their splitting, and shows the presence of the incoherent in-gap states. 
Like in Fig. 4.6, we observe zero-harmonic and higher-order modes along 
EJ + nw, which are mixed with the continuous background. In the high- 
frequency limit, the quasiparticle energies (zero-level) resemble the stationary 
Andreev states, whereas the continuum electron spectrum is far outside the 
energy gap. For such driven system, the superconducting energy gap is 
manifested in a unique manner, its properties can not be treated like in a 


static situation. 


4.4 Floquet-Magnus expansion 


In the high frequency limit, the Floquet Hamiltonian H” can be expanded 
in powers of the perturbation, called the Floquet-Magnus expansion [127,133] 


Bes yn 02 (4.57) 
v=l 
The first part is the zero component of Fourier transform of Hamiltonian 


HY” = Ay (4.58) 


0 


and the second is defined as follows 


[A + mt | fio, f,]]. (4.59) 


Here, the Fourier transform is defined as previously Hy, = + i bere H(t). 
This high frequency expansion can be adapted to the proximitized quantum 
dot Hamiltonian introduced in (3.7) with the oscillating quantum dot energy 
level 


H(t) = F (e9 + Acos(wt))dtd, + = (didi +h.c.). (4.60) 
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The Fourier transform of Hamiltonian (4.60) can be performed in few mathe- 


matical steps, presented below 
H 1 i. imut H(t) 
> "Ep e = 
T Jo 
1 T amu RA I's ot ot 
5 e (E (eh Acosta) dd + Z (dd +n) = 
1 t ia 07+ 3 1 P oe Ac, 26 3 
mi, e cadhd, +7 z |, e A; (e + e*t) dido 
LFS al een INY R ASO Taaa 
Rom, ema didi + Df, eit dnd = 
5 
00,m (Sed + sidi + Sad) + ŻT af e imwt > ę iwt(1- åt d, 
pa etm TK. 461 
zz], e (4.61) 


In the result, we obtain the Hamiltonian as follows 
tice «Danae Ween Ann? 
Hm = 50m (Send J > id) + dd) + (1m + dim) z 2) did,, (4.62) 


where the diagonal elements are connected to the discrete quantum dot 
energy level £9 and the superconducting order parameter defined by Ts/2. 
The off-diagonal elements are the result of the driving in the system. The 
first element of the Floquet-Magnus expansion, which is the zero element of 


Fourier transformed Hamiltonian, is equal to 
Fa say Ls 5 3 
HY = Hy = È eadbdo + SS dldl + Syd, (4.63) 


The block form of the Floquet Hamiltonian acquires the symmetry property 


A, = A, = 4 DA did,. The second element of the expansion can be calculated 
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as follows 
1 


mi (im Hm +e"** [ Ho, Hm] = 
m=+0 


F(2 
DE Jez > 
1 A A ROZ: A 1 X ź RAE: A 
iZ (HH + eetl Ho, H1]) + re (HH +e" Ho, Hy]) = 
„A + gimwto [ Ho, Hy) = HA, — ev imuto [Ho, H-1]) = 
[AA + eeto] Po, Ê] = eo imwto [Ho, H. 1]) = 


[AA eimuto [Ho, Hy] —e —imwto [Ho, A,]) = 


(A 
(I 
(I 
(I 


t 
hw 
t 
hw 
> (I A,,H s maton GA [Ao, H,]) f (4.64) 


This expansion requires to calculate two commutators of different Floquet 
matrix Hamiltonian (Hy, H_,] = 0, [Ho, i] = -Atsdid! + ATsdd,. Substi- 
tuting the solution of the calculated commutators into (4.64) the second term 
of the expansion is 

BO SZ [A U (45 di = AŻ 5 dt KAJ (4.65) 
In this case, the higher order terms vanish. The final expansion of the 
Hamiltonian which is truncated to the second term HE = Ho? + BE, takes 
form 


"ie Ts .ATs 
HE =Y eydyd, + (E-i Aw 


The expansion of the Hamiltonian in the high frequency regime rescales 
the induced on-dot pairing potential. The high frequency expansion results 
are consistent with the Green’s function calculations (for large frequencies) 


presented in Fig. 4.6. 


4.5 Summary 


In this chapter we studied an effective spectrum of the driven quantum 


dot placed between the superconducting and metallic reservoirs. The single 
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quantum dot energy level was periodically stimulated by the external potential. 
We have analysed the effective spectrum with respect to the amplitude A 
and frequency w of the drive. We found that the proximity induced electron 
pairing leads to the splitting of each harmonic level. The maximum distance 
between mixed states slightly decreases with larger amplitude of the drive. 
In quasiparticle spectrum, distribution of the harmonics is controlled by the 
amplitude to frequency ratio. 

The charge transport averaged over a period of driven oscillations and 
its differential conductance were also analysed. In the plots of differential 
conductance, we detected the multi-harmonic quasiparticle energies, their 
splitting and the distribution of the spectral weights. 

Additionally, we studied the influence of the finite size of the superconduct- 
ing energy gap on the quasiparticle spectrum. In a case, when the amplitude 
of oscillations exceeds the energy gap the incoherent states appear in the 
subgap regime. Those states correspond to short-time living quasiparticles 
and emerge near such values of the amplitude to frequency ratio, where the 


spectral weight of the zeroth harmonic vanishes. 


Chapter 5 


Dynamics of Andreev states in 


double quantum dot system 


This chapter presents description of the dynamical effects obtained for a 
double quantum dot system (N-DQD-S), schematically presented in Fig. 5.1. 
We consider geometry in which the quantum dots are placed in series between 
the metallic and the superconducting leads. Theoretical analysis of such 
system gives rise to the possible qubit development. The qubit based on the 
double quantum dot and similar layouts has been studied by several scientific 
groups [56-59]. Knowledge about dynamical properties of the system would 
be important for sequential manipulation and operation on them. 

In this section, we explore the time-dependent processes related to the 
creation of bound states and other properties that describe the evolution of 
the system. Specifically, we study transient effects and the evolution of the 
system after step-like or continuous changes in its parameters. In particular, 
we investigate the following quench protocols: transient effects (Sec. 5.3), 
rapid change of source-drain voltage (Sec. 5.4), quench in energy level position 


(Sec. 5.5), periodic driving (Sec. 5.6). 
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5.1 Model Hamiltonian 


The heterostructure, consisting of the two quantum dots QD; (i = 1,2) 
placed in linear configuration between the normal (N) and superconducting 


(S) leads, can be described by the following Hamiltonian 
Ê = Hs+ Hs-op, + Hpąp + HN-qp, + Hn. (5.1) 
The double quantum dot term is modelled by the single-level localized states 
Hpęb = >) Silbia + Us, + 2, (Valle + h.e.), (5.2) 
where él (ĉio) is the creation (annihilation) operator of electron at i-th quan- 
tum dot, e;, denote the energy levels, and Viz is the interdot coupling. We 


also consider electron-electron Coulomb repulsion between electrons localized 


on quantum dots U;frithi. 


QD QD, 


N 


ae — “NS 
e Dy Vio I's 


Figure 5.1: Illustration of double quantum dot's sytem placed between the normal 


metallic and the superconducting reservoirs [134]. 


The quantum dots are hybridized with the external reservoirs via 


Fin-gp, = > WTOCZ + he.) (5.3) 
ko 
and 
Hs-QD, = > (Vogt baal + h.c.) 5 (5.4) 


qo 
where Vy; (Vsą) denotes the coupling to normal (superconducting) lead. 
In absence of the metallic lead (Ty = 0) and for the superconducting 


atomic limit (Asc = oo) the Hilbert space of the proximitized DQD is spanned 
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by 16 vectors. In the occupancy representation the matrix Hamiltonian has a 
block structure, consisting of 6 subspaces [135|. Two 4-dimensional subspaces 
contain states with odd number of electrons |QD2,QD1) = |0,t), | 1,0), | t,t), 
| 1) and |0,1), | 1,0), | 14,1), | 4,14), respectively. The next two states 
| t,t), | 44) are decoupled from other ones. The remaining 6-dimensional 
subspace contains the states with even number of electrons, |0,0), |0, 14), 
1160), | 14411), | 1,4) and | 4,1), respectively. Diagonalizing the effective 
matrix Hamiltonian for U; = 0, one obtains the following eigenvalues e; and 


eigenfunctions |¢;): 


i €i lbi) 
1/2 +E a;(|0, t) F| ty, t)) + bz(=| 1,0) =| 1, 1) 
3/4 +HBFTs/2 | a,(\0,t) +] Wt) +b;(| 1,0) =| 1,14) 
5/6 +E a;(|0, 4) F|14)) +C 4,0) £| 1,14) 
7/8 +E FPs/2 a(|0,1) =| ty, 4)) + b:(| 1,0) F| J, 14) 
9 0 | t,t) 
10 0 114) 
2 (Vi2(|0,0) + | t4, 14) 
11 0 VAVŻAT$JA 
-Fs(|t,1)-|4,1))) 
12 0 all tapes) 
13/14 +T's/2 z (10,0) - | 14, 14) £10, 14) = | 14,0) 
ares Wu +I NY 
15/16 | +/4Vj5 +T'5/4 +5(10, (4) + | (4,0))+ 
Tren 14) F111) 


Table 5.1: Table of the eigenenergies and eigenfunctions of the uncorrelated 


proximitized double quantum dot system obtained for Ty =0, Ase = œ [134]. 


sa! SZ! Vi il Ei 
where E = 4 (\/4V2, +12/4+T's/2), a; = Ane and b; = 5 RTP 
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The repulsive on-dot interactions U;ń4ń,, compete with the proximity- 
induced electron pairing, leading to important changes of the bound states. 
Here, we consider their influence on dynamical phenomena within the mean- 


field approximation, using the Hartree-Fock-Bogoliubov decoupling scheme 


Ay = Prin (Rig) + Prag (Pin) + EYE (nên) + Cay Cin let et). (5.5) 
The Coulomb interactions and the superconducting proximity effect can be 
described by the effective mean-field Hamiltonian 
Hepp ® > Eig (tA Gio — X (Ai (tel él, + h.e.) + D (Viêt êz + h.c.) 
1,0 i o 


+ >> (Vurig êz + h.c.) T Y Ewkol Vp Nko, (5.6) 
ko ko 


where Eig (t) = Eio (t) = Uinio(t), A(t) = ts = Ui (ê (t)é1,(4)), 
Ao(t) = U2(ĉz (t)ć(t)). Such approximation would be reliable only in the 


weak interaction case (U <T's). 


5.2 Selected experimental realization 


The double quantum dot heterostructure has been studied intensively 
theoretically and experimentally by various groups [57, 136-138]. These in- 
vestigations established, that in this system the manipulation of spin and 
its readout are experimentally accessible [56-58]. Such possibilities provide 
opportunity to construct a qubit based on a double quantum dot, as pro- 
posed by DiVincenzo [59] for both normal leads and in Ref. [139] for the 
superconducting electrodes, respectively. 

Recently, the double quantum dot system Fig.5.1 has been further studied 
experimentally, revealing the Andreev blockade effect [64]. Mechanism of such 
Andreev blockade relies on the triplet configuration of electrons in quantum 
dots [60, 140, 141|, which suppresses the Andreev reflection processes, thus 


blocking the flow of current. This configuration offers many possibilities for 
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investigating time-dependent properties and post-quench dynamics, which 


will be addressed in this chapter. 


Al/InAs 


200 nm 


Figure 5.2: Example of the experimental realisation of the double quantum dot 


system, using InAs nanowire [64]. 


5.3 Transient effects 


We start by investigating the transient effects of the double quantum dot 
system, similarly as has been done for the junction with the single quantum 
dot [Sec. 3.4]. All results for the double quantum dot system have been 
obtained by the equation of motion technique (for details see Sec. 3.3). To 
calculate the current flowing to the normal lead jy„(t) we modify the formula 
(3.20) to following form 


tL No 
2 


Jret) = Zm Y Veet dt (êo (0)) - l, 67) 


which has been formally derived in Ref. [117]. For numerical purposes, we 
had to solve the set of coupled differential equations, which are presented in 
Appendix A.2. 

Like in Sec. 3.4, let us present the result for the evolution from the 
initial conditions to the equilibrium state. We assume that initially, at t < 0, 
the quantum dots are disconnected from both external reservoirs and next 
they are rapidly connected to the electrodes at t = 0*. Fig. 5.3 presents 
numerical results obtained for currents jy;(t), js;(t), occupancy of quantum 


dots ną(t) and the real part of the order parameter y;(t) = (d(t)dą(t)). For 
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zero source-drain voltage, we observe the evolution of the system towards 
an equilibrium state for initially empty quantum dots. In this scenario, the 
oscillations in currents and occupations are mainly caused by the flow of 
Cooper pairs into and out of the system. Those oscillations are dumped in 
time due to presence of a continuum of free electrons of the metallic lead. 
The buildup of the order parameter for individual quantum dots is presented 


in the plot of Rey;(t). Imaginary part of y;(t) is negligibly small (~ 1077). 


—Rey,(t) 


EL Teal cd ALLE 


jsą(t)[e/h] 


Figure 5.3: The transient evolution of currents jy;(t), jsą(t), occupancy ną(t) 
and the order parameter Re x;(t), imposing a sudden coupling of the quantum dot 
to both external leads at t = 0* for initial occupancy (n+(0),7,(0)). The results 
are obtained for zero source-drain voltage and system parameters: eg, = 0, Vi2 = 1, 


TN=0.2,Ts=1,U=0. 


We observe rapid development of x,(£) in time and delayed with respect to it 
formation of the order parameter in the second quantum dot. The character 


of this buildup of the order parameters directly depends on the strength to 
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the superconductor I's and the coupling between quantum dots Vj». 

Development of the bound states can by practically detected by the 
differential conductance. Using the definition of the time-dependent differen- 
tial conductance Gny(Vsa,t) = RZ JNo(t) we present the transient evolution 
of the uncorrelated system in Fig. 5.4 and the correlated system in Fig. 
5.5. These figures present the the differential conductance versus time and 
the source-drain voltage Va lifting the chemical potential of metallic lead 
(un = Vsa, Hs = 0) for various interdot couplings and the Coulomb potential. 
The transient evolution of the system exhibits emergence of the Andreev bound 
states induced by the proximity effect in the quantum dots. The Andreev 
bound states appearing in Gyo(Vsa,t) can be interpreted as the excitation 
energies between eigenstates with the even and odd number of electrons. In 
the stationary limit (t > oo) the peaks occur at E = +4 (VAV + 12/4 + Ts/2) 
(for the uncorrelated case and e,, = 0). Their broadening is caused by the 
relaxation processes on a continuous spectrum of the metallic lead. 

In the correlated case (Fig. 5.5) we observe substantial reorganization 
of the bound states. For larger U we observe smaller separation between 
their position, see panels in the first and second column of Fig. 5.5. This 
tendency proves, that the Coulomb repulsion suppresses the effects of the 
interdot hybridization and the superconducting proximity effect. In the third 
and fourth column of Fig. 5.5, the peaks merge, but presumably this is an 
artificial consequence of the mean-field approximation which is hardly reliable 


for the strongly correlated limit U > T'g. 
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Figure 5.4: The transient evolution of the differential conductance G = Gyo(Vsa, t) 
(in units of 2e*/h) with respect to source-drain voltage Vsq = LN and time t, 
obtained for Vig = 1,0.3,0.2,0.1, respectively. We used the model parameters: 
cio = 0, Py =0.1, Ig = 1, U = 0 [142]. 
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Figure 5.5: The transient evolution of the differential conductance G = Ging (Vea, t) 
(in units of 2e?/h) with respect to source-drain voltage Vsa = uy and time t, obtained 
in the weak Vj = 0.5 (upper row) and strong interdot coupling limit Vi = 2 (bottom 
row) for several values of the Coulomb potential (as indicated) for: e;o = -U/2, 


Ty = 0.1, Ts = 1 [142]. 
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5.4 Dynamics induced by abrupt source-drain voltage 


In analogy to the single quantum dot system (Sec. 3.5) we analyze dynam- 
ics caused by a rapid change of the source-drain voltage. The components of 
N-DQD-S system are initially disconnected (till t = 0) and they are coupled 


at t = 0*, exhibiting the transient evolution presented in Fig. 5.6. We observe 
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Figure 5.6: The time-dependent currents (in units of e/h) and quantum dot’s 
occupancy (a), after a sudden coupling of the quantum dots to both external leads 
at t = 0*. The post-quench evolution (b) after applying the source-drain voltage 
V;d at t = 60. Results are obtained for: ¢;, = 0, Ty = 0.2, Fs = 1, U = 0 [134]. 


the relaxation process imposed on top of the quantum oscillations of the 
characteristics. The steady state is approximately achieved at t ~ 60, after 
a sequence of quantum oscillations, damped in time with the exponential 
envelope function e”*"N/2, Frequency of the quantum oscillations coincide with 
quasiparticle energies of the in-gap bound states of N-DQD-S heterostructure. 
In particular, for e;, = 0, the period is equal to T = 47/Ts. 

At t = 60 we perform the quench, applying the source-drain voltage. 
Rapid change of the chemical potential of the metallic lead induces the 
charge flow through the system. The oscillating currents jy,(t) and jsa(t), 
are shown in Fig. 5.6b and Fig. 5.7 for different interdot coupling values. 


Here, the superconducting lead is grounded and the source-drain voltage 
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Via = UN — Hs lifts the chemical potential wy. For the source-drain voltage 
much larger than the interdot coupling |Vzg| >> Viz, we notice emergence of 
the quantum beats with the time period Tp = 7/Vi2 superimposed with higher 
frequency oscillations. The appearance of these quantum beats originates 
from a superposition of the sinusoidal waves of different frequencies. An 
analogy of this phenomenon can be found in acoustics, where an interference 
pattern results from the superposition of several sound waves of slightly 
different frequencies. In the present case for ¢;, = 0, Fy=0, U =0, the charge 
transport involves the Andreev scattering related to position of the in-gap 
bound states, which are formed at energies E = +3(\/4V,3 +T2/4+Tg/2). 
It has been shown [120-122], that for the single quantum dot coupled to 
the normal leads, the jump of the source-drain voltage induces coherent 
oscillations in the current characteristics, with frequency w = |Vsq — €g|. For 
the double quantum dot setup, we should replace the quantum dot energy 
level by the four quasiparticle energies at which the Andreev scattering is 
amplified. Taking into account the number of these quasiparticle energies, the 
total current is a superposition of sinusoidal waves X4; a;e~t sin(Q,t), with 
frequencies (ją = Veg £ wą and Q3/4 = Veg + wą, Where wyją = Vio + I's/4. The 
sinusoidal waves are damped with corresponding parameters A; and multiplied 
by the coefficients a; which control the the contributions to total current from 
these in-gap bound states. 

In Ref. [143], the authors have shown that in the normal junction the 


beating pattern depends on the ratio 


_ Wy, + W? (5.8) 


© [wy = wl 
which can be applied to our case. This ratio can be expressed by the coupling 


between the dots and hybridization to superconducting lead 


AVi> 
=——, 5.9 
rao (5.9) 


The numerical results of jy,(t) obtained for V,a = 20 (Fig. 5.6b) show that 
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Figure 5.7: The post-quench currents jyo(t) (in units of e/h), obtained for several 
values of the interdot coupling Vj2 and source-drain voltage V:a, using: Eio = 0, 


Ty = 0.2, Ts =1, U = 0 [134]. 


the ratio is r = 8 and the following sequence of the beats can be detected: 
7,7 x 5. If the ratio r is not an integer number, then the beating patterns 
are significantly more complex. The plot of jy.(t) current shown in Fig. 
5.6b indicates, that for large Vq = 20 the post-quench evolution is mainly 
composed of the quantum oscillation with period Tp = 7/V\2 superimposed 
with the faster oscillations, whose frequency is equal to Vsa. In the steady 
state limit, the current jy„(t) is larger for Vsa = 2 than for Vq = 1.5, because 
of the broader transport window which involves all in-gap Andreev bound 
states. Additionally, we notice that the current js,(t) differs from jyo(t), and 
their difference is more pronounced upon increasing Vsa. We can attribute 


this behaviour to the fact that quantum dots between the external leads wash 
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out small fluctuations of the current js,(t). It also imposes the period of 
oscillations equal 4r/T's. A closer look on this beating structure presented in 
Fig. 5.7, displays the time-dependent current jyo(t) after the abrupt rise of 
the source-drain voltage. Here, we clearly see the influence of the interdot 
coupling Viz on the interference pattern by observing the oscillations with time 
period Tg = r/V;>. In real measurements the period of beating oscillations 
could serve as a tool to evaluate the interdot coupling Viz. For example, if the 
coupling to the superconducting lead T's ~ 200ueV and Viz = 0.5I's, I's and 
21's the beating period would be Tp ~ 21, 10 and 5 picoseconds, respectively. 


This time-scale is nowadays attainable experimentally. 
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Figure 5.8: The same as in Fig. 5.7 but for the correlated half-filled quantum dots 
U=0.5. 


Fig. 5.8 displays results obtained for the correlated case assuming both 


quantum dots to be half-filled £; = -U/2. In this particular case, the evolution 
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after the quench looks almost the same as in Fig. 5.7. We notice only the 
slightly different shapes of the curves. Nevertheless, the oscillation period 
remains unchanged. The different profiles of the time evolution between 
uncorrelated and correlated cases can be explained by almost the same initial 
conditions right before the jump of the source-drain voltage. In both cases, 
the functions of quantum statistical averaging of operators (...) have almost 
the same values for half-filled scenario. In other cases than the half-filled 
scheme, we can expect analogous behaviour of the oscillations if initially 
(before the jump) the quantum statistical averaging of observables takes the 


same values with and without correlations. 


5.5 Quench in the quantum dot energy levels 


In this section, we explore dynamics of the double quantum dot system 
driven by the change of the quantum dot energy levels. We consider the 
sequence of quenches imposed on the quantum dot energy levels ¢;,. First 
quench occurs at tı = 60, when the energy level of both quantum dots are 
raised up by the gate potential Vj: Eio > Eio + Vg. Second quench takes place, 
when the system achieves a new stationary state (at time ty = 120), then we 
abruptly change the energy levels back to their initial values: e;, + V} > Eo. 

Let us first inspect the quantum dot occupancy of the stationary limit 
t > oo (Fig. 5.9) of the quantum dots with respect to the initial energy 
level €1, = €9,, for a few values of the interdot coupling Viz. We notice that 
quantum dot occupancies are almost identical nj, ~ n2,. Here, we observe 
step-like functions of the charge occupancy characterized by plateau regions. 
Upon varying the energy level the charge evolves from the filled to nearly 
empty quantum dots, through the plateau regions corresponding to ny, z 1, 
0.5, and 0, respectively. The width of the half-filling plateau region depends 
on the interdot coupling and we observe its contraction with smaller Vj». 


Changeover between plateau regions occurs at Eis x =V1o. 
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Figure 5.9: Charge occupancy of QD; in the stationary limit as a function of the 
energy level €2, = €1, determined for several interdot couplings V12. The dashed line 
refers to the correlated system, U = 1. Results are obtained for: V,g = 0, Ty = 0.1, 
Ts = 1 [134]. 


In Fig. 5.10, we present the current jyo(t) obtained for the strong interdot 
coupling Viz = 4 and identical energy levels e;, = 0. We display the post-quench 
evolution after the first and the second quench for several gate potentials 
V, = 3.0, 3.2, 5, and 6, respectively. We notice that the post-quench evolution 
has the same general properties if the energy level position e,, corresponds to 
the same value of n;,(t = oo) for arbitrary Viz. We thus need only to consider 
the case of any value of the interdot coupling and a set of various jumps of 
Eigi 

In Fig. 5.11, we show na,(t), jno(t), Jso(t) obtained for the interdot 
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Figure 5.10: Evolution of jya(t) (in units of e/h), after step-like variation of the 
quantum dots energy levels e; > Eio + Vg at t = 60 and Ej, + Vg > Eic at t = 120. We 
assumed model parameters: €;, = 0, Vi2 = 4, Vesa = 0, Ty = 0.2, Ts = 1, U = 0 [134]. 


coupling Vız = 4 and energy levels ¢;, = 0, where the post-quench evolution 
after the first and the second quench had been obtained for gate potentials 
V, = 3.2, 3.8, 4, and 5, respectively. Such values of V, correspond to the 
quantum dot’s occupancies equal to ~ 0.48, ~ 0.4, ~ 0.25 and ~ 0.015 (see 
Fig. 5.9). In the first case, the jump from e; = 0 to Eio + V} = 3.2 (top rows 
in Fig. 5.11), we notice that the stationary value n,(t ~ 120) = 0.48 and it 
is nearly the same as the initial quantum dot occupancy n2,(t = 60) = 0.5. 
Therefore, m2,(t) remains nearly unchanged exhibiting only small fluctuations. 
Likewise, the currents jy, and js, exhibit negligibly small oscillations in 
their characteristics. After the second quench, the occupancy ngs is almost 
intact ~ 0.5, but immediately after the quench the transient beating structure 
appears. For the larger amplitude V, (V, = 3.8 and V; = 4), the fluctuations 
are more pronounced. We observe the oscillations with period T = 7/V\2 and 
the beating structure with the period Tp = 2r/T s. Rising the value of gate 
voltage V} up to V} = 5 and above, the evolution after the quench differs from 
the previous cases. After the first quench, na, (for V, = 5) it shows almost 


exponential decay, down to nearly zero value with the damped oscillations. 
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Figure 5.11: Evolution of the quantum dot’s occupancy (a) and currents (b-c) (in 
units of e/h), after step-like variation of the quantum dots energy levels Eo + Eio + Vg 
at t = 60 and e;, + Vg > Eig at t = 120. We assumed model parameters: €;, = 0, 


Vio = 4, Veg = 0, Py = 0.2, Is = 1, U = 0 [134]. 


In the evolution, after the second quench, the oscillations of quantities with 
the period T = 2r/T s are visible without the beating structure. We conclude 
that, after the second quench, the beating structure appears if e,, corresponds 
to the slopes of nay(t = 00) with respect to £2 (Fig. 5.9) and it is absent if 
Eio coincides to the bottom plateau of these same curves. For Viz < 1, the 
evolution after the first quench is similar to the case with large Vj2, but after 
the second quench the beating patterns are not visible, only the oscillations 
with the period 47/I’s can be detected. 

We have also performed calculations considering the Coulomb interactions 
U =1 and applying the same quench procedure, see Fig. 5.12. For V} = 3.35, 
4.0, 4.25 and 5.3 (which correspond to the tilted part of the curve calculated 
for U = 1 in Fig. 5.9), the occupancies ne, are the same as in the uncorrelated 
case. In the post-quench evolution we obtain almost identical behaviour of all 
quantities. Like in previous chapter (Sec. 5.4), the initial conditions before 
the quench are almost the same for uncorrelated and correlated cases, which 
is the reason why the post-quench evolution is almost identical. This specific 


scenario presents the weak influence of the correlations on the post-quench 
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Figure 5.12: The same as in Fig. 5.11 but for the correlated system U = 1, assuming 
the initial energy levels ¢;, = -U /2. 


evolution. Notice that the results in Fig. 5.12 were obtained for U =I's5 on 
the borderline of applicability of the mean-field approximation. For larger 
Coulomb potential the calculations need to be done by more sophisticated 


methods. 


5.6 Periodic driving of orbital level positions in DQD 


We extended our investigation of driven systems, discussing the double 
quantum dot system under the drive of the energy levels ¢;,(t) = Asin(wt). 
For simplicity, we assume that both quantum dots oscillate with the same 
amplitude A and frequency w. Left panel in Fig. 5.13 presents the time- 
dependent current js,(t) obtained for w = 0.1 and for several amplitudes A, 
neglecting the correlation effects. The right panel in Fig. 5.13 presents the 
time-dependent current js,(t) for various Coulomb potentials U. Top panel 
of Fig. 5.13a shows the transient evolution of the unbiased system (Vsq = 0) 
which is characterized by damped oscillations with the time period 47/T's. 
Similar behaviour of the time-dependent currents is visible in the next plots 
presented in the left panel. We observe different profiles of the quantum 
oscillations damped in time (up to t = 20) and the long-term fluctuations 


caused by the periodic driving appearing upon increasing amplitude A. In the 
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Figure 5.13: The time-dependent currents (in units of e/h) for the periodically 
driven e;„(t), where the oscillations are turned on at t = 0*. The panel a) presents 
the results obtained for Viz = 4, U = 0 and several amplitudes A. The panel b) 
presents the results obtained within the mean-field approximation for Viz = 3, A=3 
and several values of U. The dashed line shows the profile of the oscillations energy 
levels (not in scale). The other model parameters: Va = 0, w = 0.1, Ty=0.1, 
Ts=1 [134]. 


case of strong interdot coupling, Viz = 4, we notice that the time-dependent 
currents vanish in the asymptotic limit for amplitudes A < 3.5 (see Fig. 5.13a). 
For A > 3.5, the amplitude A exceeds the energy of subgap quasiparticles 
therefore the oscillation of the time-dependent currents can be observed. It 
turns out that the current js,(t) remains oscillating through all time region, 
including the asymptotic limit t > oo. We observe that the period of this 
oscillatory behaviour matches with the driving period. 

In the right panel in Fig. 5.13, we present the numerical results for the 
current js,(t) obtained for the Coulomb correlations U; = U2 = U within the 
mean-field approximation. For comparison, the upper curve in Fig. 5.13b 
presents the results for uncorrelated case. The results are obtained for the 
unbiased system Vq = 0 with interdot coupling Vi2 = 3 and the driving 
amplitude A = 3. We note that this case, Viz = A, corresponds to the situation 


when the current js (t) oscillates in the entire time window. 
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Figure 5.14: The time-dependent currents (in units of e/h) for the periodically 
driven e;„(t), where the oscillations are turned on at t = 0*. The results are obtained 
for the different periodic driving schemes with profiles marked by the dashed lines 
(not in scale). The other model parameters: A = 2, w = 0.1, Vig = 4, Vea = 0, 
Ty = 0.2, Ts=1, U =0 [134]. 


Fig. 5.13b displays the influence of electron correlations. With increasing 
U, we observe the suppressed amplitude of the current. However, the details 
of the oscillations remain still visible. 

Studying the periodically perturbed double quantum dot heterostructure 
we now consider different profiles of the drive. Fig. 5.14 presents the time- 
dependent current js,(t) imposed by the periodically driven energy levels 
Eig(t) = €i¢(t+ T). In each plot, we included the corresponding driving profile 
marked by dashed lines (the amplitude is not in scale). Moreover, we assumed 
the same amplitudes and frequencies. For the step-like variations of the 
quantum dot’s energy level, we observe the enhancement of the current after 
each change of e,,. Those oscillations are damped in time and their lifetime 
is smaller than the period of the drive. For more smooth profiles of ¢;, the 
amplitude of the revival oscillations decreases and the current characteristics 


are correlated with the current plot for sinusoidal changes of the quantum 
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dot’s energy level. 

In the periodically driven system, we can calculate the averaged quantities 
over one driving period, e.g., current and differential conductance. This has 
been done for one single quantum dot system in Sec. 4.3 using the Green’s 
function method. The observable averaged over one driving period gives more 
precise information about the role of the amplitude A and frequency w of the 
oscillating quantum dot’s energy level. Here, we use the Heisenberg equation 
of motion method, inspecting the charge currents averaged over a single period 
T = 2m/w of the perturbation and we analyse the quasiparticle features visible 
in the nonequilibrium transport properties of the investigated heterostructure. 
We calculate the Andreev current average over one driving period (jno(t))t, = 
A a jNa(t)dt and its differential conductance Gyo(Vsa) = a ino (t))to 
with respect to the source-drain voltage V,a and driving amplitude A or 
frequency w. 

Till t < 0 we assume both quantum dots to be empty and (at t > 0) impose 
the driving on e;, simultaneously with applying the voltage Vig = Un — ps. 
The time-dependent evolution of the currents is presented in Fig. 5.13, where 
we observe the oscillations due to the driving force. In a long time evolution, 
the current becomes nearly periodic in time, therefore we can calculate its 
average value. It is important to notice that the period of the drive and the 
period of the currents become the same in the long time limit. In Fig. 5.15 we 
present the numerical results of the differential conductance Gyo(Vsa,w) for 
the set of the model parameters, where ly «I's and ca = 0. The maps present 
the averaged Andreev conductance obtained for two values of the interdot 
coupling Vj. and several amplitudes A, as indicated. Panels in Fig. 5.15 
display the characteristic features of the driven quantum system originating 
from the absorption/emission of the external field quanta. In general, we 
notice the main quasiparticle peaks appearing at +3 (/AVŻ + 12/4 +T s/2) 
and their higher harmonic peaks. Moreover, the avoided crossing behaviour 


between the peaks of different harmonic levels can be noticed. 
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Figure 5.15: The averaged over one driving period Andreev conductance Gyo (in 
units of 2e?/h) as a function of frequency w and the source-drain voltage V;g. The 
results obtained for several amplitudes A and interdot coupling Vig (as indicated). 


The other model parameters: Iy = 0.1, Ts = 1, U = 0 [134]. 


Some properties of the photon-assisted tunnelling through the quantum 
dots placed between the normal electrodes have been previously studied in 
literature [144-146]. For such system, the main resonance peaks and harmonic 
bands are modulated by the square of the Bessel functions of the first kind 
J?(A/w) (see Fig. 4.3). In the Josephson junction, where the single quantum 
dot is embedded between two superconducting leads, the quasiparticle peaks 
are proportional to the following Bessel function J?(2A/w) [69,147]. Here, 
the argument is doubled due to the fact that the Cooper pairs are the charge 
current carriers. 

In our N-DQD-S heterostructure, we observe that the harmonic peaks 


are weighted by the squared Bessel function J?(2A/w). For the constant 
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Vo=1.0, w=1.0 V = 2.0, w=1.0 
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Figure 5.16: Gy, (in units of 2e?/h) averaged over one driving period Andreev 
conductance as a function of the amplitude A and the source-drain voltage Vsa- 
The results are obtained for several frequencies w and interdot coupling Vj2 (as 


indicated), using: Ty = 0.1, Is = 1, U = 0 [134]. 


amplitude A = 4 (see Fig. 5.15d), the quasiparticle peaks disappear at 
frequencies w for which the Bessel function is equal to zero, e.g. w ~ 3.3, 
1.45, 0.92. They correspond to the first, second and third zero of Jo(2A/w), 
respectively. Moreover, comparing Figs.5.15c and 5.15d, we see that the value 
of w for which the main quasiparticle peaks and higher harmonics disappear 
is independent of the interdot coupling V2. 

To investigate the influence of the amplitude on the quasiparticle spectrum, 
we calculated the averaged Andreev conductance with respect to V,a and A, 
for a few values of the interdot coupling Viz and frequency w (see Fig. 5.16). 
In the stationary situation, for A = 0, there exist four peaks in the differential 


conductance. Upon increasing the driving power, the main quasiparticle peaks 
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loose part of the spectral weights (intensity) at expense of the higher order 
replicas. With the growing amplitude, the larger number of replicas can be 
detected at energies w, +w, +2w. Simultaneously, their total spectral weight 
undergoes substantial redistribution. Besides this tendency, for specific values 
of A, the spectral weight of individual harmonics vanishes and then again 


reappears. 
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Figure 5.17: The averaged Andreev conductance Gyo (in units of 2e?/h) as a 
function of the interdot coupling Vig and the source-drain voltage Vsa. The results 
are obtained for several amplitudes A (as indicated), using the model parameters: 


w=1,Ty =0.1,Tg =1, U = 0 [134]. 


Fig. 5.17 presents the plots of averaged Andreev conductance against 
the source-drain voltage Vsa and the interdot coupling Vi», obtained for two 
values of the driving amplitude A and frequency w = 1. The left panel in 
Fig. 5.17 presents the quasiparticle peaks appearing around +nw. Upon 
increasing the coupling Viz the peaks gradually split into the lower and upper 
branches. The branches of conductance peaks never cross each other because 
of the quantum mechanical interference [148]. In the case with the larger 
driving power (right panel in Fig. 5.17), we also observe this avoided crossing 
behaviour. Moreover, the harmonics consist of two nearby located peaks. 

Finally, we investigate the influence of the electron correlation on the 
differential conductance Gy,. Fig. 5.18 shows the results for Gy, as a 


function of the source-drain voltage V,, calculated for several values of the 
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Coulomb repulsion force U and the following set of parameters: Vio = 2, 
A = 2 and w = 2.5. In comparison to the first plot presented in Fig. 5.18 
(for uncorrelated case), we observe two doubled main quasiparticle peaks 
localized at V.j = +E* = +5 (AVG + 12/4 + Ts/2). We observe the first-order 
harmonics at V,a = +E* F hw and the peaks from the second-order harmonics 
localized at V,qg =+E* F2hw. For larger electron correlation, the quasiparticle 
peaks change their heights. Upon increasing U, we observe more asymmetric 
distribution among the second-order harmonics as compared to the first-order 
harmonics, which have nearly symmetric shape. Position of the quasiparticle 


peaks are nearly insusceptible to the electron correlations. 


Figure 5.18: The averaged Andreev conductance G = Gy, (in units of 2e?/h) as a 
function of the source-drain voltage Vsa. The results obtained within mean-field 
approximation for several values of U (as indicated), using: Vją = 2, w = 2.5, 


Ty = 0.1, Ts = 1 [134]. 
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5.7 Summary 


In this chapter, we discussed the transient effects and the post-quench 
dynamics of the double quantum dot system. We presented that the transient 
effects are characterized by damped quantum oscillations with the exponential 
envelope function e~'/?, The period of these quantum oscillation coincides 
with energies of the in-gap states for the system. 

The post-quench evolution induced by jump of the source-drain voltage has 
shown that the quantum oscillations are the superposition of the sinusoidal 
waves 4, ae sin(Q;t), with adequate frequencies Qij2 = Vsa +w and 
Q3ją = Vsa + wą, Where wyją = Vio + Fs/4. For the large source-drain voltage 
[Vsa] >> Via, the beating pattern depends on the ratio r = 4V;2/T's. 

We also analysed evolution of the system after the sudden shift of the 
quantum dot energy levels. We noticed, that the evolution of the system after 
the quench mostly depends on the initial occupations n;, (right before the 
quench) and the magnitude of gate potential V}. For the energy shifts which 
coincide with the plateau and slope regions presented in Fig. 5.9 (plot of the 
charge vs. energy level of quantum dots) some universal properties of the 
oscillations can be detected. 

In the last part we presented the results for the system influenced by the 
periodic driving of the quantum dot energy levels. The oscillations of the time- 
dependent currents are induced when the amplitude of the drive is close or 
exceeds the value of the coupling Viz. We noticed that for stronger correlations 
the amplitude of the oscillations decreases, but the characteristic shape of 
currents remains unchanged. We also presented the conductance averaged 
over one driving period. Such conductance is characterized by the main 
quasiparticle peaks appearing at +5 (AVG + 12/4 +Ts/ 2) and emergence of 
the higher harmonic peaks caused by the emission/absorption of the energy 
quanta. The correlations affect the distribution of the quasiparticle peaks, 
mainly by their height (spectral weight). Position of the quasiparticle peaks 


weakly affected by the electron correlations. 


Chapter 6 
Machine learning simulations 


In this chapter we discuss a computational method of a machine learning 
(ML) [149,150] applied in our study to investigation of the driven systems. 
Considering the double quantum dot system, the numerical solution of the 
Heisenberg equation of motion consumes a large amount of computational 
power and time. Our numerical procedure has been parallelized which greatly 
increased the effectiveness of our code. However, to obtain a single map of the 
averaged conductance over one driving period using CPU 2x Xeon E5-2660 
(2.2GHz 16 cores/ 32 threads) took almost one week of calculations. To avoid 
such difficulties an auxiliary procedure based on a machine learning algorithm 
has been developed. 

In this section the concept, main principles and structure of the machine 
learning algorithm will be introduced. For brevity we present here only the 


most important details. 


6.1 Fundamentals of artificial intelligence 


In general, the machine learning is a numerical method of data analysis. 
The purpose of the machine learning is to expose patterns and correlations 
between the particular variables of the investigated case. The idea of the 


machine learning algorithm has been presented in 1959 by Arthur Samuel [149], 
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who was the pioneer in an artificial intelligence (AI) field. The machine 
learning is a subfield of the artificial intelligence, which can be regarded as a 
simulation of the human mind. The main purpose of the artificial intelligence 
is to mimic humans actions, and program machines to behave like humans. 
Those simulations can be applied to learn and solve some specific problems. 

The artificial intelligence is based on an artificial neural network, which 
is composed of the single nodes called perceptron (see Fig. 6.1). Moreover, 
the idea of preceptron is based on the working principles of the neuron. In 
general, the preceptron can have any number of inputs and outputs, which are 


represented by numerical values. Every input is multiplied by a corresponding 


Figure 6.1: Idea of the perceptron. 


weight. Then the result of the multiplication is passed to an activation 
function and finally, in the output, the specific result is obtained. Sometimes, 
the activation function requires additional bias value. It happens when the 
sum of inputs does not fulfil the properties of the activation function. For 
example: The perceptron returns 1 if the sum of inputs is positive number or 
0 if the sum of inputs is negative number. The problem occurs if the original 
input is zero. Then an addition of the extra bias value equal +1 solves the 
problem. Mathematically, the concept of the perceptron can be pressented 
by a simple sum yy wizi + b, where N denotes number of inputs, w; is the 
weight corresponding to the specific input z, and b is the bias value. 

A larger number of perceptrons can create the matrix form that is the 
neural network, schematically presented in Fig. 6.2. The particular parts of 


the neural network can be specified: input layer, hidden layer, output layer. 
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The input layer consists of the perceptrons whose number usually corresponds 


input hidden output 
layer layers layer 
_— 


Figure 6.2: Graph of the artificial network. 


to the number of features or properties of the considered case. In the center 
of the neural network, the hidden layers are located. The number of the 
hidden layers and the number of nodes in the specific layer is an arbitrary 
value. Often the hidden layers are called the level of abstraction because 
their outputs do not represent any physical value. At the end of the neural 
network, the output layer is located. The final layer returns the predicted 
value of the formed model. 

In the model, besides the number of layers and nodes, a crucial role plays 
the activation function. Depending on the task, the appropriate activation 
function must be applied. If we consider a classification problem then the 
logic function will be the appropriate choice. For more complex cases the 
more dynamic functions will be better. The most common functions having 
good performance are: Sigmoid function - S(x) = (1+ e-*)-!, Hyperbolic 
Tangent - tanh(x), Rectified Linear Unit (ReLU) f(x) = z* = max(0, z). 


6.2 Basics and functions of the learning protocols 


To grasp the idea of the learning process, some particular information 
needs to be presented. First of all, let us consider the case of the single 


perceptron and its predictions. 
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For random weights, the output of the perceptron and the expected 
value are inconsistent. Adjustment of the perceptron weights is called the 
learning process. To measure the discrepancy between predicted value and 
the expected value we can use a cost function. The cost function C evaluates 
the performance of the neuron that is measured by how far off the predictions 
x are from the true value y. One of the most common function is a quadratic 
cost C = Yav(y: - 1;)2/N, which is not an efficient function and could increase 
the time of the learning process. For this cost function, the larger errors are 
more pronounced due to the squaring. Another, more efficient cost function is 
a cross entropy C = -4 Yo [yi n(z;) + (1- y,) ln(1- z;)]. The cross entropy 
allows for a much faster learning process because the larger the difference, the 
faster the perceptron can learn. The cost functions are indirectly dependent 
on the perceptron weights. 

A subsequent important aspect of the learning process is a gradient descent 
(the delta rule). The gradient descent is an algorithm for finding the minimum 


of the cost function C, schematically presented in Fig. 6.3. As we see in 


C(w) 


Gradient 
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Figure 6.3: Visualization of the gradient descent algorithm. 


Fig. 6.3 there exists a value of weight w which for the cost function C(w) 
has its minimum. The gradient descent algorithm is based on an iteration 
process. The algorithm starts with a random set of perceptron weights and 
calculates the error for each next step of the iteration, updating the model 


parameters until the result is the minimum of the cost function. This process 
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has following pseudocode structure 


while: C(w,) #min[C(w)] { 
Ow; i 


Wi = Wi-a 


where a denotes a learning rate. The learning rate define how quickly the 
gradient descent algorithm moves towards the minimum of the cost function. 
Inappropriate value of the learning rate negatively affects the gradient descent 
algorithm, if a is too small the algorithm makes small steps that extend 
calculation time. On the other hand, if a is too large the algorithm can 
overshoot and do not find the solution. One can specify three types of the 


gradient descent algorithm: 


e Batch gradient descent - all of the training data is considered by an 
algorithm to make a single step. The algorithm takes the average of the 
gradients of all the training points and then uses that mean gradient to 
update the perceptron weights. One step of gradient descent is called 


one epoch. 


e Mini-batch gradient descent - the algorithm takes small portion of the 
training data called batch. The batch is used to update the perceptron 


weights in each iteration. 


e Stochastic gradient descent - updates the perceptron weights in each 
iteration using only a randomly selected single batch from the training 


data set. 


In more complex cases, when the neural network is considered, the larger 
number of perceptrons and their weights are needed to figure out the solution. 
To adjust the values for the perceptron weights and minimize the cost function 
across our entire network by the gradient descent algorithm we need to use 
a backpropagation. The backpropagation algorithm is a generalization of 


the gradient descent algorithm for the artificial neural network. Such an 
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algorithm calculates the gradient of the error function with respect to the 
neural network’s weights. The calculations proceed backwards through the 
neural network. At first, the algorithm calculates the gradient of the loss 
function for the final layer, then going back calculates the gradient of the 
loss function using the weights of the subsequent layers up to the first layer. 
The backpropagation algorithm is an algorithm for a supervised learning. 
Supervised learning is one of the types of machine learning, it uses the labeled 
data sets to model dependencies between the target predictions and the input 
features. 

The learning protocol is a combination of the particular procedures pre- 


sented above, schematically visualized in Fig. 6.4. To train the neural network 


Figure 6.4: Scheme of the learning protocol for the artificial neural network. 


we require a prepared dataset, which will be divided into three subsets: train, 
test, holdout. The train subset contains the data for the training process, 
which teaches our neural network by using the gradient descent and backprop- 
agation algorithms. Then the efficiency of the established neural network is 
checked by data from the test subset. If the predictions of the model are not 
convergent to the true values the training procedure will be performed again. 


The final correctness of the neural network is verified by the holdout subset. 
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6.3 Adaptation of artificial intelligence to the driven 


double quantum dot system 


In our specific case of N-DQD-S heterostructure, we train the neural 
network using supervised learning and the mini-batch gradient descent. Our 
data set is a compound of 76 conductance maps. The maps of the conductance 
averaged over single driving period have different resolutions and in total give 
971760 data points. To increase the number of points, every map has been 
linearly interpolated, giving us 3887040 data points in total. As software 
we choose Tensorflow, the open-source machine learning library having an 
integrated programming interface Keras [151]. 

In the model of the densely connected neural network (see Fig. 6.5), we 
have four parameters in the intake layer (Viz, w, Vsa, A) and the output layer 
with the single neuron (G), where the hidden layers are composed of the 


layers consisted of 2048, 1024, 512, and 256 neurons, respectively. The size 


hidden layers 
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Figure 6.5: The densely connected artificial neural network. 


of the neural network is caused by the nonlinearity in the system. Every 
single neuron has the sigmoid function as the activation function. To avoid 
overfitting of the model we defined a dropout of 1% neurons on every learning 


step. In the training process, we established the batch size 1024 and the 
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number of epochs equal 600, which gave us the fidelity coefficient R? = 0.987. 
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Figure 6.6: a) The plot of the differential conductance predicted by the neural 
network versus calculated values. b) The reconstruction of the conductance map in 


Fig. 5.16a generated by the neural network [134]. 
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Figure 6.7: The calculated conductance map (panel a) and generated by the neural 


network (panel b) obtained for Vig = 1.7, w = 2.5 [134]. 


v. 


Fig. 6.6a shows the discrepancies between the calculated and predicted 
values of G for our model, which is deflation of the points from the red 
line. Fig. 6.6b presents a reconstruction of the conductance map for specific 
system parameters, which has been included in the train data subset. For 
comparison, panels in Fig. 6.7 show the conductance maps obtained from 


the direct calculation (panel a) and by the neural network (panel b). The 
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map shown in panel a) in Fig. 6.7 has not been used for learning the neural 
network. Those results have been published in [134]. 

In Fig. 6.7, we see the qualitative reconstruction of the calculated by 
equation of motion method conductance map. With this approach, we are able 
to simulate new conductance maps for a periodically driven double quantum 
dot system and many different values of the controlling parameters. Recalling 
that we can control: source-drain bias voltage V4, amplitude A and frequency 
w of driving field and the strength of interdot coupling Viz. The analysis and 
properties of these maps are discussed in Chapter 5.6. The machine learning 
model of the double quantum dot system allows us to simulate the same 
conductance map faster. We have managed to reduce the calculation time to 
less than an hour, where calculations using the equation of motion method 
required more than a week. 

This neural network model is available at the following webpage: 


www. dropbox.com/sh/Ohzs9im3d3bf0 jr/AADRr3kltw2m0dCCh8tedolWa?d1=0 


Chapter 7 


Summary and Outlook 


This dissertation presented a study of the superconducting nanostructures and 
their time-dependent properties. We have investigated the systems composed of 
the single and double quantum dots embedded between the normal metallic and 
superconducting leads. We have inspected response of such systems on various 
dynamical perturbations, i.e.: transient effects (sudden coupling of the quantum 
dots to the reservoirs), abrupt application of the source-drain voltage, quench in the 
quantum dot energy level, sudden change of the coupling to the superconductive 
electrode. Moreover, we have analysed the evolution of the system under the 
periodic driving of the energy level of the quantum dots. Variety of the protocols 
has resulted in a plethora of emergent physical phenomena. Here, we briefly recollect 


the most important findings. 


Single quantum dot 
For the single quantum dot system, placed between the normal metal and 
superconducting electrode, we observed the time-dependent emergence of the An- 


dreev states. Formation of the bound states in the quantum dot at energies 


E4 = +4/ (eq + U/2)? + (Ts/2)? has an oscillating character, with period 2r/E4 and 


the relaxation processes are governed by envelope function e^t. 


After the rapid switching the source-drain bias voltage in the uncorrelated 
quantum dot system, we also observed the oscillations of the charge current. Those 


oscillations resulted from the superposition of two contributing currents with the 
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frequencies w, ją = LN +T 5/2. We noticed that the current to normal lead jNo(t) is 
consisted of the high frequency component w = u and the low frequency one with 
the time period T = 2r/T s, whereas the current jsa(t) revealed oscillations with 
only the lower frequency. 

In the post-quench evolution of uncorrelated system (imposed by variation of 
the quantum dot energy), the oscillatory behaviour of the observables revealed 
the period T = 2r/ li + 4(Aeg)?. For the correlated quantum dot system, we 
discovered the phenomenon resembling the 0-7 transition, sign-reversal of the 
oscillating observables. Investigating the influence of the suddenly changed coupling 
to the superconducting electrode, we observed the transition from the doublet to 
singlet configurations. 

For the uncorrelated system we found that the periodical driving field induced 
the harmonic states in the quasiparticle spectrum. Such harmonic states emerged 
from the Andreev states and showed up in the differential conductivity of N-QD-S 
heterostructure. We also noticed that the harmonic states linearly depend on the 


driving frequency and they split due to the superconducting proximity effect. 


Double quantum dot 

For the double quantum dot system placed in series between the normal 
metal and superconducting electrode, we observed the time-dependent forma- 
tion of the Andreev states appearing in the differential conductance at E = 
+2 ( /AVĄ +T3/4+ rs/2) (for e; = 0). 

In the post-quench evolution induced by the source-drain voltage, we noticed 
emergence of the quantum beats with the time period Tp = 7/Vi2. The appearance 
of the quantum beats is due to the superposition of waves of different frequen- 
cies. The total current can be represented as a superposition of sinusoidal waves 
> 3 aye" sin(Q;t), with frequencies Oj = Vsa + wą and (03), = Vsa + wą, where 
Wij2 = Vi2 Ts/4. For the correlated case at the electron-hole symmetry point 
Eig = —U/2, we observed similar behaviour of the observables. 

Subsequently, we investigated the quenches imposed on the quantum dot energy 
levels €;,. We established the protocol, consisting of the first (ei, > Eio + Vg) and 


the second quench ( Eio + Vg + Ego) for Eio = 0. After the first quench, for V} < Vio, 
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we observed the damped oscillations in the characteristics of observables and almost 
unchanged occupancies of both QDs. For V} ~ Vi2, we observed more pronounced 
oscillations with the period T = 2m/Ts and variation of the occupancy of quantum 
dots. In the evolution, after the second quench, for V} ~ Viz, the fluctuations were 
more pronounced. We observed the oscillations with period T = 7/Vi2 and the 
beating structure with the period Tp = 2r/T s. For V} > Viz, the evolution after 
the second quench revealed the oscillations with the period T = 2m/Ts without the 
beating structure. The post-quench dynamics proved to be strictly dependent on 
the ratio between V} and V12. For the correlated case, at half-filling Eis = -U/2, we 
observed similar behaviour of the observables. 

Additionally, for the double quantum dot system with the periodic driving, we 
found appearance of the harmonic peaks weighted by the squared Bessel function 
Ję(2A/w). They originated from the Andreev bound states and evolved depend- 
ing on the amplitude and the frequency of the driving field. These branches of 
conductance peaks never cross each other. For investigating the emergence of the 


harmonic peaks we constructed the machine learning model (discussed in Ch. 6). 


Outlook 

In the future studies, we would like to investigate Majorana modes in systems 
with the periodic driving. In our opinion, the Majorana modes should have a 
impact on the conductance map of the system with the single or double quantum 
dot. In multi-terminal systems with the side-attached Majorana wire |152, 153], 
the dynamics induced by the different types of quench protocols would be worth 
addressing as well. Additionally, it would be interesting to explore the influence of 
strong electron correlations on the post-quench dynamics of states in the double 
quantum dot. 

Furthermore, it would be worthwhile to investigate the dynamics of quantum 
dot arrays connected to superconducting reservoirs, and analyse their evolution 
induced by different types of perturbations. Another issue for investigation could be 
the time-dependent entanglement between quasiparticles in heterostructures with 


the topological superconductors. 


Appendix A 


A.1 Differential equations of motion for N-QD-S 


Solving the Heisenberg equation of motion for N-QD-S Hamiltonian we get the 


closed set of equations: 


di fi = de(nay(t)) = 21m p Ve ENE fer — sp R A. (A.1) 
k 
di fo = di(na,(t)} = 2Im b Ve "NE! fea — iS p = A. , (A.2) 
k 
di f3 = dild, (t)dy(t)) = (e(t) + Eq (t) = ily) fs = iAsc(1 m f1 = fo) 
-i J. Ve NE fra + iY Vee NE fies, (A.3) 
k k 


dyfka = dr(di(t)cyę(0)) = GO = | Jka + Ase fka 


+iVpe Nt (Ek), (A.4) 
di fun = dy(d (#)emy(0)) = [ica (t) $ x| TeRi 

+iVpe ENE (ex), (A.5) 
thea ala [iea Š | EOT (A.6) 
difia = delay (Eer (O) = [iea (t) = | far + isofin, (A.7) 


where n(ex) = [1 + exp((ex - UN)/kBT)] " is the Fermi distribution function of the 


normal lead electrons and d; = 4. 
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A.2 Differential equations of motion for N-DQD-S 
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The closed set of the differential equations corresponding to the N-DQD-S 


system Hamiltonian: 
di fi = dz(ny(t)) = 21m [Vio fz — Ascfs], 
dt fo = dz(na, (t)) = 2Im [Via fg — Asc fs], 
di fz = dz(naę(t)) ) = 2| 5 e` Ekt fyg — EN p- Viel. 


Vk 
difa = dę(na,(t)) = 21m 2, Vie" ‘eel fig — ZE Vial. 


di fs = di(dy,(t)diy(t)) = — len) + ea, (t)] fs + iVi2 [ fo — fio] 
-iAsc[1- fı- f2], 


di fe = dz(doj (t)day(t)) = ~i [ean (t) + ex(t) — iI N] fe + iVi2 [fo — 


+i > Vee "FH" [fer - frs], 
k 


fio] 


(A.12) 


(A.13) 


dy fz = dy(di, (t)dzy(t)) = IO = ex(t) + w fr + iVi (fs = fi) + iAscfio 


é -igit 
-i X Ve ** fir, 
k 


(A.14) 


di fg = dild} (t)d2(t)) = IO -Ex (t) + ra] fs + IVq> (fa — fo) — iAscfo 


-i Y Vie fho, 
k 


(A.15) 


dę fo = di(diy(t)da,(t)) = —i [eno + €2, (t) — w fo + Wa [fs + fe] — iAscfs 


-i Ve firs, 
k 


(A.16) 


di fio = dz(du(t)dy(t)) = —i MO + ezę(£) — wy fio — iVi [fs + fe] + tAscf7 


: -igit 
-i X Vie ** frs, 
k 


(A.17) 
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di fra = deldi, (t)crt(0)) = ieqy(£) fer + Vio fis + Ase fea, A.18 
di fr = dild? i (t)cxy(0)) = ten (t) fka + Viz fro — Asc hrs, A.19 
dt fra = dz(dy(t)cz,(0)) = —terp (t) fea — Vi2 Ser — iAsefk2, A.20 
di fra = de(diy(t) cx (0)) = -tery(t) fra — Vio tas + Ase fia, A.21 
di frs = ded, (Ł)cjy(0)) = 


+iVpe""'A(Ex), A.23 
di fke = dę(dź, (t)cr,(0)) =2 BO + rx] fke + iVi2fk2 A.24 
+iVie" i(ex), A.25 


(A.18) 
(A.19) 
(A.20) 
(A.21) 
ilro + fra] fks + Vio fki (A.22) 
(A.23) 
(A.24) 
(A.25) 
(A.26) 


difur = de(day(t)c44(0)) = -i|ezr(t) = SP | far = ‘Vio, 
di fas = di(da (Pe (0)) = ~i | e210) - ra] Sst. AI 


Such systems of differential equations can be solved by iterative methods. In 
our investigation, we applied Runge-Kutta 4th order method. However, other and 


faster numerical methods, e.g. Dormand-Prince, can be implemented. 


A.3 Non-equilibrium Green’s function formalism 


In quantum physics the many-body system can be governed by time dependent 
Hamiltonian H (t). The idea of adiabatic switching relies on assumption that in 
the past t = —co the system was in initial state defined by the many-body density 
matrix p(-00) and from that point evolves according to the Heisenberg equation of 
motion d p(t) = -i[A(t), p(t)]. For h= 1 the formal solution is represented by the 
initial state and the evolution operators p(t) = U(t,-00)p(-0o)Uf(t, -00), where 
the operators are defined as follows 

O(t,t’) E Tel H(r)dr _ Jim 6 1HBÓt,-IH(t-6t)6t o-i (t +6t) dt (A.28) 
and ôt = (t-t’)/N denotes the infinitesimal time step. The expectation value of 
observable © for the system at specific time t is defined as 
Tr{Op(t)} _ THO (00,1) OU (t, -00)p(-»0)) 


(OO) = GO} TO) | 


(A.29) 
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where the trace is performed over the Hilbert space of many-body Hamiltonian 
[154-156]. The numerator of the fraction includes the time evolution operator 
describing the system evolution from t = -o0 to time t and backward to t = —-oo. 

The zero temperature expectation value of the operator O in equilibrium is 
defined as (GS|O|GS) = (0|U(-oo,t)ÓU(t,-00)|0), where |GS) = U(t,-00)|0) is a 
ground state of the interacting system and |0} denotes a vacuum |157]. The evolution 
operators describe the evolution of the non-interacting vacuum state |0) toward the 
ground state |GS) through the adiabatic activation of interactions. It is claimed 
that the adiabatic evolution of the system along the entire contour and back to 
the starting point results in an additional phase factor (0|U (+00, —00) = (O|e’”. To 
compensate the phase factor in the definition of the expectation value for operator 
we can divide by this term (GS|O|GS) = (0|U(+00,t)OU(t, —00)|0)/e”, then the 
evolution describes only the forward path without backward segment. 

This approach does not work in nonequilibrium situations. For the case, when 
the system was driven out of equilibrium, the final state depends on the switching 
protocol and the entire evolution of the system. The final state does not have to be 
the same as the initial state. To avoid the needed of knowing state of the system at 
t= +00, we can use the Schwinger approach [158]. The evolution of (A.29) can be 
extended from t to t = co and backward by substituting U(t, +00)U(+oo,t) = 1 and 


written as 


(Ó(Ł)) = eee (A.30) 


which is the evolution along the closed time contour C. In temperature limit T > 0, 


those two branches of the time evolution are depicted by the Keldysh contour in 


Fig. A.1. 
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oo 


ton Cx tN+1 


Figure A.1: The Keldysh contour C including the forward C+ and the backward C- 


branches. 


In the Keldysh formalism, where the time evolution of observables is governed 
by the contour, the path-ordered Green's function is defined as a quantum statistical 


averaging of the fermion operators 
GGC, t) = -i(Tofai(t); (t’)}), (A.31) 


where ae ) is creation operator for electron in state j and time t’, ć;(t) is annihila- 
tion operator for electron in state į and time t. The single particle Green’s function 
describes its propagation through time and space. Here To is a time ordering 
operator which governs the operator’s positions on the Keldysh contour [154]. The 
contour-ordered Green’s function can be divided on four different types, depending 


on the forward or reverse branches: 


i To {éi(t);e1(t')}) tł eC, "casual", 


Gylt) = -i( 
G (t,t) = -i(E;();E!(Y)) tEC_,t' EC, "greater", 
GE (t,t") = (tt) = -i(ć;(t);ć,(£ )) + "g 
GY (t,t) = +i(E!(E);ćz(t)) teC,,t'eC_ "esser", 
-i( 


GS (t,t) = T.(ć;(t); é(t')}) tt eC_  "antitime-ordered". 
(A.32) 
The "casual" Green's function is also called "time-ordered" Green’s function. For 


the fermionic operators the time-ordering operator Tc fulfils the following conditions 


A SHË) ift>t 
To{e&(t); E= E (A.33) 


AANA: W, 
Ae a(t) Et >. 
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Fermionic operators obey the anticommutation relations: 


{ĉ;, ĉj} = 0, (A.34) 
600, (A.35) 
12,0) = Oi, (A.36) 


where anti-commutator is defined as 44, B} = AB + BA and 6;; is the Kronecker 
delta. 

Having defined the contour-ordered functions, we can write down fundamental 
objects of the Keldysh technique: the retarded (r), advanced (a) and Keldysh (k) 


Green’s functions 


G;, (t,t') = -i6(t — t')({é(t); E(E'))) 
=P(ft-t)[G;(tt') -CKE t], (A.37) 


Gę(t,t') = 10C — t)(46;(t); ie) 
=0(t - t) [C t) - Gź,(t,t)], (A.38) 


Gź(t,t') = GY) + GGE), (A.39) 


where 8(x) is a Heaviside step function. These relations are only valid when t # t’. 
Since by definition of the greater and lesser Green’s function in (A.32) [G<]! = -G*, 


one notices that 


[a"]! = G°, (A.40) 


[GF]! =-G*. (A.41) 


For non-zero temperatures T > 0, the Green’s function becomes more complex. 
In order to obtain information about the particle and its space-time position as 
a function of temperature, we need to introduce the Matsubara Green’s function 
[159, 160]. Thermal averaging (Ó) = Tr (e22) = Tr (ee?) = Tr (tir) is 
equivalent to evolution along imaginary axis from 0 to -i(, see Fig. A.2. Therefore, 


the Matsubara’s approach is called the imaginary time approach with 


t > it. (A.42) 
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t=t forward 


Figure A.2: The Keldysh or Konstantinov-Perel’s contour in the complex time plane 


with forward C,, backward C_ and imaginary branches. 

Recalling that, for general operators A and B, the Green’s function is homogeneous 

with respect to time, i.e. it depends only on a difference of the imaginary times: 
GAB(T1,72) = Gap(t1 — 72) = GąB(T), (A.43) 


and that the function G4g(7) is periodic (anti-periodic) function for bosonic 


(fermionic) Green's functions 
GaB(T — m$) = -Žr (P fA (r- m3) B) (A.44) 
with a general argument 7 being in range 
mB<r<(m+1)8, (A.45) 


where m denotes an arbitrary integer, we can show its other properties. If 7—m > 0 


the Green’s function fulfills 
1 ik 
Gap(t - mB) = Gir (eP" BA(r - (m+ 1)8)), (A.46) 


and if 7 < 0 the operators obey BA(T) = —nTc¢ A(r)B, where 7 is an integer number. 


Since T - (m + 1)5 <0 the Green’s function follows 
Gą4B(T - mB) = —qGaB(T-(m+1)5). (A.47) 
In particular case, when m = —1, we see that 


GuB(T) = -7Gap(T + 8) (A.48) 
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and the Matsubara functions are periodic (anti-periodic) for bosonic (fermionic) 
choice of the time ordering. Since both functions are periodic with period 28 the 
Matsubara Green’s function can be expanded in the Fourier series as follows 
In > ; 
GĄpB(T)=g 2, e” GiBlwm) (A.49) 
m=- o0 


with frequencies distinct for bosonic and fermionic propagators 


a for bosons, 
Üm = (A.50) 
Cm Dz for fermions. 


A.4 Langreth rules 


In the uncorrelated case the contour ordered Green’s function satisfies the Dyson 


equation of the general form 


Gt) =g Et) f dr | drg TETT), (A.51) 


where g°(t,t') is bare Green’s function and X°(r,7') a self-energy term. In this 


equation the time convolution between the functions occurs via 
C(t,t') = i dr A(t, 7) B(r,t'), (A.52) 


with the integrations over the extra time argument 7 on the Keldysh contour. We 
can assume that t is on the first part of the contour and t’ belongs to the second 
half. The time convolution between the functions deforms the Keldysh contour. For 


future investigation let us analyse a "lesser" function 


C“ (t,t) = Ja drA(t,r)B'*(r,t') + | dr AT (t;T)B(qt'). (A.53) 


Lesser function C< (t,t) can be split into two integrals with two different branches 
Cı and C2, which means that we deform the Keldysh contour by adding two 
extra branches. In the first term, the sign < for function B occurs because in the 
integration over contour C; argument 7 is always lesser than t’ (similar for the 


second term). The first part of (A.53) can be slit into two parts 
t —00 
f, drAtnB<Gt)= fo dra (7) B (rt) + | | drA'(t,r)B<(rt') 
1 —00 t 


= f ” drA'(t,r)B<(r,t') (A.54) 
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due to definition of the retarded Green’s function (A.37). The second term on the 
right hand side of (A.53) can be analysed in the similar manner and in the result 


B function is "advanced" function 
I drA'(t,r)B(r,t') = 1 drA'(t,r)Be(r,t'). (A.55) 
2 2: 


Combining (A.54) and (A.55) we can rewrite (A.53) in form of the first Langreth’s 


rule: 
cx(tt)= f dr[A"(t,1)B (r,t) + A(t) Br], (A.56) 


The "retarded" or "advanced" component of a contour ordered C (t,t) (A.52) 


function is 


O(t, t") = J: is dr A(t, 7) B”? (r, t"). (A.57) 
The proof for "retarded" function is as follows: 
C" (t,t) = 8(t-t)(C7(t,t')-C'*(t,t')) 
- 0-1) [| dr[A" (t, TAB? (r,t) - B<(r,t')) 
+ (4° (t, 7) - A*(t,r))B"(r,t')] 
-60-5[[_dGlPl7)- AGB (r,t) - BC t) 
+ f dr(A? (t,r) - A(t, (Bt) - B*(r,t'))] 
å f, arar(e Br t’). (A.58) 


For the time convolution of more than two functions C(t, t) = fo ABC the "lesser" 


or "greater" function type of C(t,t') posses the form 
O? = T [A" B'C>< + A" B><C* + A><BOC"], (A.59) 
where the integration is over real time axis. The lesser function C< reduces to 
one J. A’ BSC, (A.60) 


when a steady state is turned on by adiabatic evolution from t = -o0 to t. 
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A.5 Convolution in Fourier space 


With the Keldysh contour function C(t,t’) can be rewritten in the form of 
product of other functions C(t, t’) = [% drA(t,r)B(r,t'). Below, the proof of the 
time convolution between two-times function in Fourier space is presented: 


c(t,ty= | drA(t,r)B(n,t) = | dr F 


kymy 


w/2 . ; 
i: deze ere eerie) A n (€1)x 
RE -i(E2+lw)T+i(E2+moaw)t" 
`y J deye ‘*? POU Bima (E2) = 


l,ma w/ 


w/2 s 3 ; . r 
( Il 5 [ is] eo tlerthw)tti(ertmiw)7 „-i(e2+lw)Tr+i(E2+maw)t Akm (€1) Bim. (e2) z 


i=1,2 Ni, Mj 


w/2 s 4 r 
| II b> a is] gle thaliri(earmat)t 0(€1 — E2 + M1 W — lw) Apm, (€1) Bim, (€2). 


i=1,2 n;m; 2 79/2 
(A.61) 
Then we can make the following substitution: 
le € (-w/2;w/2) > ó(e1 —€92 + (m = Dw) = 0(€1 = €1) Om, I] 
obtaining: 
w/2 } ; A 
Cat F J , dee i(e*ko)tri(etmaw)t 4, (e) Bim, (€). (A.62) 
k Lm Y ~#/2 
In the next step, we change the summation index as follows 
|k > n,l > k, ma > mı| to obtain 
w/2 ; NE Y 
C(t)= ¥ A . dee Keer A. (2) Bupske)s (A.63) 
nym,k ~~ 


After the Fourier transformation, C(t, t’) becomes Cym(e) equal to the convolution 


of two functions in Fourier space: 


Cnm(e) = Y Ank(€) Bem (€). (A.64) 
k 


A.6 Bare Green’s function in Floquet space 


By definition, the bare Green’s function for quantum dot has following form 


Gi, at, (ot) = -ilo {dolt d(t)}). (A.65) 
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Solution of the equation of motion for the bare Green’s function 9, dł (t,t’) of the 


single quantum dot described by Hamiltonian H(t) = X„ea(t)dyd, requires the 


time derivative of operator de(t) 


Adz(t) = iL H(t), do(t)] = ie o WHO) H(t), doje WHO) (A.66) 
= jej lo WEGO A(t), dz Je | fo WEW) (A.67) 
= R edgr (tet fo WHO) [at (t)dor(t),do(t)Je to WEW) (A.68) 
= jeg: (thet fo WHO) g et ho WEW) _ ie 4(t)d,(t). (A.69) 


In result, we obtain 


tig’ p (tt!) = be(t-t' boo - (To {idide(t), di (t')) = (A.70) 
Selt- t' boo" — t€a(t)(Tof{do(t), dt, (t’)) = (A.71) 
Se(t—t' bao! +Ea(t)g, p (tt), (A.72) 


which can be recast to form: 
(i, calt)) 9, p (tt!) = delt- door (A.73) 


where €q(t) = £9 + Acos(wt) is time dependent discrete energy level of quantum dot, 
e9 is the stationary quantum dot energy level and A amplitude of the perturbation. 
Writing the discrete time dependent quantum dot energy level eg(t) = Èk gee eM 
in Fourier space and applying the Langreth’s rule (A.57) we obtain the retarded 


and advanced equation of motion form of the Green’s function 
(ið; + i0* — >LETNa (t,t) = Selt- t')5o,9". (A.74) 
k o> o! 


Here, our aim is to obtain time-independent equations. For this purpose, we apply 
the Fourier transform to both sides of the equation 


: w/2 ; ; ' 
(ia pie >| T [ 7 dee tetnw)ttilermw)t gr (e) = 
k nym 77W 


da,di „mm 


>] a a ae (A.75) 
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Rewriting the equation in the form 


w/2 l ? , 
(id; re i0*) > da dee *Etw)tri(ermv)t goe (e)- 


da,di „mm 


w/2 , 
3 = dee '(ete)t+i(etmw)t o -iwkt SA g’* (e) = 


s 
Amk do,d_,,nm 
» R dee '(Etw)tri(etmo)t 5 
Z Jij nm) 

and substituting the following calculated derivative iett) = (c+ nw)e 


we obtain: 


w/2 iat -i(c+nw)t+i(e+mvw)t' r,a 
2 a de (e + 10” +nw)e g (e)- 


dg dt „nm 


w/2 : ; 
- k tak , 
) : ie dze i[e+(n+k)w]t+i(e+mw)t CI at mE) = 


3 ios dee tetnw)tti(ermu)t! z 
n,m -w/2 ~_ 


We introduce the extra summation over the additional index k 


3 Ta de (e Eat + kw) eo tletkw)t+i(etmw)t! Ón.k kd i p. (e)- 


n,m,k 


4 
-ife+(k+k)w]t+i(e+mw)t! k ra A 
>: Jin? Sanh at, am lE) 


> ead dee *(Etw)tri(etmo)t z 
woja nm) 


and replace the index k by k =n - k to finally get 


n,m,k 


Edo 9 


w/2 
-i(etkw)t+i(e+mw)t! en k ra = 
1 dee dod}, 3,2 z 


n,m,k 


DEE i deg here re 2 


> pe de (e + i0* + kw) ę tletkw)t+i(e+mw)t! Önk M d, = (e)- 


132 


(A.76) 


-i(c+nw)t 


(A.77) 


(A.78) 


(A.79) 


(A.80) 
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In the next step we can skip the Fourier integration on both sides and write the 


result in the form of time-independent equation 


OE ue one en” |g” a „dł, mom! 


> 


r,a-1 


Sde? = (e = i0* + nw)ópm — Ego”; (A.81) 


which is the bare Green's function in the Fourier space. 


A.7  Diagonalization of bare Green’s function 


In the preceding part, we obtained the energy-dependent Green's function in 


Fourier space (A.81). Here we introduce A function [129], defined as follows 

Kom = (lonktom)) = (ep |. dteMetu(t)) = lola f adtelr-™"eltto(t)) = 
L | dów mech (o|U(t,0)|o) = AE dtci(n—m)ut „iegt e-i fy dreda(7) — 

> T jenot o-i ją dr(edo(1)-c (A.82) 


where ((on| acts on a state in the Floquet space |uo,.m)). The operator A fulfils the 


following property 

Aa nA! NE » i x |, ater Bate dr(edo(r)-=9) 1 = | demek dr(Edo(T)-E3) — 
dz z T dy, 8 _i fo 

ie ln- k)aę + fo dz(edo (z/w)-€9) (a i(m- k)yę i Sy dz(edo (2/w)- e9) — 

fs dx i as: eilne- my)-ik(x- ve- ó Jo dz(Edo(z/w)- ed BU fy dz(edo (z/w)- e9) = 


+ On 2 


fs da % dy ele E ie dz(€qq(z/w)-e9 = 
> k 


J.= dx NIE ról - yje- my) rere dz(Edo (z/w)-€9) = 

T na. = Oni: (A.83) 
-r 2T 

Here, we used the mathematical property of the Dirac delta f(x) = f ô(x-y)f (y)dy. 


The product of two A operators gives the Kronecker delta $% Ne gh = Onm. To 
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obtain the diagonal kind of the bare Green’s function D' Ą d (e) we act on both 
do, o nn 
ra-1 


di, kl 


pr 1 r,a—1l 
A dt m AOR Dank „di, MOUE = 


FA a (e + 10* + kw) x1 — Edo.k-1]Aaln = 
kl 


sides of the inverse Green’s function 9, (e) by A operator 


2 AV gle + i0" )Sk Aon + 2 Ni kwon No.in = » NATE LA = 
2 At ACE 10°) No kn + 2 AL „KOŃ okn = 2 AM eye Ney 
E + 10° + (nw + Ego )Onn — Edo,n-n — Edo,k-IOklOln = 
e i0" + nw — Edo,k-lÓkl = € + 10" + nw - Ede, (A.84) 
where we used the following relation 
> donkkwh! pm = (nw + e9) See (A.85) 
k 


To prove (A.85), we need to rewrite the limits of integration in (A.85) in the 


following manner 


Di Manik pm = [73 = [oy aS lna my)-ik(x-y) bije ees dz(Eqc(z/w)-€9 


pe dx =f dye i(ng- MY) kwe” F ET eg), (A.86) 
Using the definition of the Dirac delta and its differential 


1 -ik(xz- 
d(x-y)= > ie ie? 
k 


10,0(a — y) = i(-1)kd(a - y) =kó(xz-y), (A.87) 
we perform some analytical calculations, to obtain 
Aner Te Zf dyei("®-™9) ins] O,5(a - y)]e 24 laelle) -D — 
k 
enrm) 55 (x = yee Jy dz(Edo(2/w)-eg) pes 


J. =f dye GM) S(x — y) [nw - eao (z/w) + Eye" hr zlEdsElw)-cg) = 


i dE ilma yy - Edo (z/w) + eg] = * dE i(n-m)e( 
-r 20 -r 2T 

i dZ i(n-m)e 
-r 20 


nw + £9]- 


Edo (z/w) = (nw + £9) — £do,m-n- (A.88) 
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In result, the inverse diagonal bare Green’s function is given by 


ra-l a z: 
D dt nn © = € i0" +nw- €go,0; (A.89) 


where nw — €go,9 is the eigenvalue of Hamiltonian Hp in state |usn)). Then the 
inverse diagonal Green’s function can be written as 
D; (e) =e 10° - ((uon|Hrluan)). (A.90) 


dod, nn 


In a similar way, we can obtain the bare Green's function expressed by A operators 


and determine the diagonal part of the Green's function in the Floquet space 
„> „dł oO A onkD dz p publ) Ao, kniz eT) 


For the oscillating energy level eq(t) = £9 + Acos(wt) the operator Agnm takes the 


form 


Ao,nm z 


a dte i(n-m)wt eile dr(eg+A cos(wr)-e9) _ 
1 i(n-m)wt ,-i fg dr(e0+Acos(wT)-s0) 1 t i(n-m)wt „—i JĄ drAcos(wT) 
= f° ate e "Jo OT" DD = f dte e "Jo OT pa 
0 0 
OCZ JAM dh: 
> ję dti motei sin(wt) (A.92) 
where we used the Bessel function property [161] 


ec iasin(wt) = » eet] (a) (A.93) 
k=-00 


Ao, NE ave dte' (= m)wt > P (=)- 
A 1 ; A 
a dt i(n-m-k)wt y (>) Rz TS dt AES i(k-k)wt y (=) = 
T A > A Va) TR > WE A 
1 fT = A A 
> dt Ôk n-mJk | — | = Jn-m | =]. A.94 
T J 2 : 5 ( w ( w l ) 


to express the operator A by the Bessel function dependent on the Fourier coefficients, 


the driving amplitude, and the frequency of the perturbation 


A 
Agnm = Jn-m (=) . (A.95) 
w 
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Now we are able to transform the bare Green’s function as follows 


(e) = EAD KRONE 


a, nah „um 


F Joan (5), do „dj ew ©) Te- „(5) o 


Jn-k (5) Irm (s) 


, A.96 
3 E+i0* + kw- E9 a) 
and get the final formula 
oo = ; A. 
dh gn) 25. E+ i0° + kw- 5 (E 


A.8 Mixed Green’s functions 


Quantum dot - metal interface 
The contour-ordered Green’s functions which describe the quantum dot coupled to 


the metallic lead are defined by 


Ge at, (tt!) = -i(fotcto(t) $), (A.98) 


Cko, 


GĄ a, (ht) = {Tole O EG). (A.99) 


To determine these Green’s functions we must solve the following differential 


equations 


i0,GQ at (6t) = -ilfo {kcen lt), dh, (t)}), (A.100) 


iG | w (tt!) = -ilfo {acho (t),dL(t)}). (A.101) 
Nko’ o! 


To obtain the solution of such functions we need to calculate the equation of motions 
for the operators CNko(t) and cho (t) with the corresponding Hamiltonian H(t) 
defined in (3.1) 


1O:CNko (t) = il H(t), CNko(t)] = ENko(t)CNko(t) + Vkdo(t), (A.102) 


iðichpo (t) = iLH (E), chro (t)] = -Ento (tcro (t) —Vedl(t), (A-103) 
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which give us 


(iði — ENko(t)) Gai, Mt te VG at, (© t"), (A.104) 
(iði + ENko(t)) G „i, bt) = -VeGa at (6t). (A.105) 


On the left side of (A.104) and (A.105) there appear the bare Green’s functions 
of the normal lead electrons defined by 


9, t (t,t') = -i(Tof{enko(t), Cypor(t’)}), (A.106) 


CNko Cup gl 


gy (ht) = -ifeleke lt), enro (YP. (A.107) 
Nko'ĆNko 
For the normal metallic lead Hamiltonian A N = Dko € NkoC\ ne C Nko the bare Green’s 


function obeys 


(i0;-Enko(t))go + (t,t) = 6-(t-t) doe", (A.108) 
CNko npg! 

(iði + Enzo (t)) 9% (t,t) = felt- t') bao". (A.109) 
CNko:"Nko'! 


Substituting (A.108) and (A.109) to Eqs.(A.104),(A.105), respectively, we obtain 


dad (t, r) 7 Vi Je WARE (t, T)G at (7, t’), (A.110) 
c ACZ i x : 
Ga, CE) = Va [ATI (b DG, at Tt) (A-111) 


Quantum dot - superconductor interface 
Below, we present the equation of motion for the contour-ordered Green’s functions 


for the quantum dot coupled to the superconducting lead, described by Hamiltonian 


Hs = ie ee = Żyj Ag C + h.c.) 


(iô: m Esq) E t) = AG dł (t, t) F VąC at (t, t’), (A.112) 
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. c PY: c iS OSI € 1 
(iô: T Esq) RAZ ) = ARE as ) Va at di he ). (A.113) 


This set of entangled functions comes from the solution of the equation of motion 


for operators 


idicsqo(t) = EsqCSqo(t) + Ojo Ascch q(t) - Ste Asch q (t) + Vadot),  (A.114) 


ich qo (t) = -Esqch_za(£) — So AscCsgt (t) + Sto Ascesa (t) — Vedi (t).  (A.115) 


To solve functions presented in (A.112) and (A.113) we need a solution of the 


following bare Green’s functions 


‘Deg at (tt) = õelt =t) - i(fotidcsgo(t), Czgo:(1'))) (A.116) 
CSq sah 
idege a (tt!) = elt = t’) - i(TofiOesqo(t), qoi(t)}), (A.117) 
CSąt:CSgt 
which give 
(ið; = Esq) gF i (t,t') = õc(t z t’) = Asci i (tt), (A.118) 
CSqatCSqt “s-qlSat 
GO + Esa) a (t,t )=-Ascg® + (tt). (A.119) 
CS-a Sat CSatCSqt 


Similarly, we can obtain a set of two equation for g° +, (t,t). In such way, we 
CSalCSqy 
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can obtain other required bare Green's functions 


(a)? - EZ) g a (t,t) = Ascóe(t-t'), (A.120) 
CS-qt'"Sql 


((i0,)” TAD t (t,t) z -Ascóe(t - t"), (A.121) 
CS-a Sat 


where Eg, = Eka FA. 
Expressions for the bare Green’s functions, (A.112) and (A.113) can be trans- 


formed as follow 


GS t,t’) =-V, [ drg® 
pee ) the ŚLĄ 


c c 1 
+ J: dro, EVG (rt) (A.122) 


c I 
ba (t, T)G at (7, t ) 
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G° t,t =-V; | drg LDCS (r,t 
ch „at , ) q c Tica , ) hatl ) 


c c 1 
-V, Jos, ogy IG g(t ). (A.123) 


A.9 Contour ordered self energies 


This appendix presents calculations of the self energies for electrodes. To solve 
the Green’s function presented in (4.39) we need the self-energy functions in the 
Fourier representation. 

We start from the self-energy of the normal metallic electrode 


(7,7') = AA A nr), (A.124) 


a Ionro c 


which can be Fourier transformed in (4.32) into the following form 


xi O E Ve (A.125) 


CNo Cy, nm k CNko:;CĆNko 1 EE = + nw — Ek 


where ék = ENko — EF and ep is the Fermi level of the lead. By changing the 


summation over the momentum vector into integration we obtain 


oo V? pp. 
2d gps ft Ce (A.126) 
CNo Cg NM —00 e+40* + nw -— Ep + in 


where px is density of states of the electron in the metallic lead. Using the 
relationship for 7 > 0 


| 1 
lim > 3 
n>0e+10* + nw — ék + in 


= p.v. Finð(e + i0" + nw- &), (A.127) 
where p.v. is the principal value, we get in the large bandwidth limit 
5%, (e) = Finónm A GEV 58 (E240 Paw = Ee: (A.128) 
CNo;CyG: NM — 00 
Assuming constant coupling for every momentum vector k: Vk = Vy and pk =p, we 


obtain 


>. (€) = Fin VE pónm. (A.129) 


CNo;CHy „um 
In the wide-band limit approximation, we obtain the constant coupling defined by 


T = 2rpV?, which finally gives us 


T 
Dii „(e) = F bam: (A.130) 


T 
CNo CNG 
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To find the self energies of the superconducting electrode, we start from 


er = AŻ a (ar), (A.131) 
S-qt' Sq 


"a 


where the bare Green’s function has the form 
((60:)°- E3,)9% + (tt) = Ascde(t- t’). (A.132) 
S-qt? Sal 
From the double Fourier transform of (A.131) we obtain 


Ae) = Vasi (e), (A.133) 


teh nm ch ar bap nm 


B 
with 
Asc 
r,a scónm 
; €)= 3 A.134 
Ś ean ) (e +107 + nw)? e ść = AŚ l ) 


By changing the summation over the momentum vector into integration, we obtain 


PqAscônm 
d = 
ô 
= 2A scp V2 d ża 
BRB Ją EED + nw)? —-@- AD, 


A nm [O(Age — |e + nwl) F isign(e + nw)O(|e + nw| - Ase)] 
yl(e + i0* + nw)? - AŻ.| 


= Vénp 
(A.135) 


In the wide-band limit approximation we obtain the off-diagonal self-energy 


yra Ce PsAsc Onm [O(Ase — lE + nw|) F isign(e + nw)O(|e + nw| - As.) 


cbpebpnm 2 VIE # i0* + nw)? — A? 


(A.136) 


In similar way, we can obtain the diagonal self-energy of (4.41). We use the Fourier 


transformation of g” A t") 
Sach 


(£) = V9 s (e), (A.137) 


im CSqsC tm 


with 
ô 
o e) = L = A.138 
R ) (e + 10* + nw)? =: © = A2 ( ) 
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Again changing the summation over the momentum vector into integration, we get 


00 + nw)ó 
ag -00 Sa Va (e = i0* + nw)? - £2 - A2, 


me ô 
- oV2 Í (e + nw)ónm 
PVs 0 da (e +10* + nw)? - ¿2 - A2, 
nm [O(Asc — |E + nw) F isign (e + nw)O(|e + nw| — Asc)] 
Lr 
a OG (Cees 


(A.139) 


In the wide-band limit approximation the diagonal part of self-energy takes the 

following form 

yra oe -Tg(e +nw) drm [O(Ase — |e + nw) F isign(e + nw)O(le + nw] - Ase) ] 
i = 


CSC nm 2 yl(e + i0* + nw)? - A2] 


(A.140) 
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2. Transient effects in a double quantum dot sandwiched laterally between super- 
conducting and metallic leads, R. Taranko, K. Wrześniewski, B. Baran , I. 


Weymann, I. Domański, PHys. REv. B 103, 165430, 2021 
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